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PRECEDli^Gi 


SUMMARY 

The primary objective of this study was to provide a data base for 
a program plan- for the development of the ion-propulsion thrust system 
for the Halley's comet mission spacecraft. This data base was to include: 
the definition of a design concept, selected from among alternate candi- 
date configurations; the identification of required supporting technology, 
including the definition of critical areas and potential technical risks; 
the definition of a program development plan, including a development 
■schedule and an assessment of potential schedule risks; and a preliminary 

estimate of yearly and total program costs . 

A concurrent objective of the study was to conduct a hardware 
"approach confirmation" technology effort to evaluate the ion thruster s 
performance and lifetime at the power level required for the Halley's 
comet mission, to design and evaluate the thruster isolator required for 
operation at the higher power level, and to evaluate the design of a 
capacitor-diode voltage multiplier. 

A thrust system baseline configuration was identified for the 
30-cm extended-performance mercury ion thruster that can perform the 
Halley's comet rendezvous mission. The configuration is comprised of 
10 thrusters configured with a power management and control system and 
, a structure and thermal control system in a modular thrust system design. 
The power management and control system uses conventional power process- 
ing. Power is provided to the thrust system with an 85 kW concentrating 
solar array. The thrust system mass is 1010 kg (including 15/4 contin- 
gency), the average system efficiency is 70%, and the estimated relia- 
bility upper bound is 72%. 

Adaptability of the 900-series 30-cm thruster design to the 6 to 
7 kW range required for the Halley's comet mission was demonstrated with 
only minor design modification required, and an acceptable high-voltage 
isolator design was validated by laboratory tests. The design and per- 
formance of an alternate power management and control system design 
approach utilizing the capacitor-diode voltage multiplier was successfully 
demonstrated by laboratory model tests at power levels in excess of 1 kVJ. 
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The technology efforts mentioned above assisted in the identification 
of the level of technical risks associated with the thrust system design. 
These risks have been found amenable to resolution through normal engi- 
neering development and, therefore, judged to be acceptable for mission 
application. 

The program plan, which includes the procurement plan generated for 
the baseline configuration, is a viable plan that provides for delivery 
in May 1981 of the flight thrust system to be integrated with the mission 
module and solar array. The cost of the thrust system development pro- 
gram is projected to be 54 million dollars (in fiscal year 1977 dollars) 
excluding contractor fee, of which approximately 13.5 million dollars 
will be required in fiscal year 1978. 

In contrast to the low technical risk, the schedule risk for 
initiating this program development is of particular concern. Timely 
approval of the authorization of 13.5 million dollars for fiscal year 
1978 must be granted so that the pre-project, or advanced development, 
activities can be initiated. 
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SECTION 1 


INTRODUCTION 

This report summarizes* the results of a six-month study to define 
the design, program plan, and costs of the ion-propulsion thrust system 
for the Halley's comet mission spacecraft. The modular characteristics 
of the design developed during this study also make it applicable as the 
prime space propulsion system for other potential missions. 

This study, which is based on an initial system characterization 
(completed 7 February 1977) performed by the National Aeronautics and 
Space Administration's Lewis Research Center (NASA LeRC), was performed 
in three parts: 

9 Design tradeoff studies (14 February to 15 April 1977) 
to define and compare alternate design approaches 

• Conceptual design definition, program plan, and costs 
of a selected design approach (15 April to 15 June 1977) 

6 Approach confirmation of supporting technology in 
selected areas. 

The results of this study are presented in five volumes.. Volume I 
summarizes the results of the entire program. Volume II discusses the 
conceptual design, program development plan, and cost estimates for the 
selected baseline thrust system design. Volume III describes the design 
tradeoff studies- performed to compare alternate design approaches. 

This volume. Volume IV, describes the evaluation of thruster technology 
for extended-performance applications. Volume V presents the details 
of the capacitor-diode voltage multi pier (CDVM) circuit analysis and 
experimental evaluation. The results reported in these volumes have also 
been presented in briefings at NASA LeRC. 

A. BACKGROUND 

In the fall of 1976, the Office of Aeronautics and Space Technology 
(OAST) was given the responsibility of assessing the capability of the 
electric propulsion technology under development at NASA LeRC and of the 
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solar array technology under development at Marshall Space Flight Center 
(MSFC) and the Jet Propulsion Laboratory (JPL) to perform the Halley's 
comet rendezvous mission proposed by JPL. OAST established an "August 
Project" team from members of the three organizations to develop a 
preliminary program plan to support a fiscal year (FY) 1979 new start. 

The August Project consisted of parallel efforts by JPL, NASA LeRC, 
and MSFC to define the design approach, program plan, costs, and risks 
of the Halley's comet mission. Three areas were considered: the space- 

craft (including the science payload), the ion propulsion subsystem 
(referred to as the thrust system in this report), and the solar array. 

The NASA LeRC program was conducted in two phases. First, initialization 
studies (completed 15 February 1977) were conducted to define require- 
ments and to identify preliminary design characteristics. Second, during 
the 15 February to 15 July period, alternative thrust-system configurations 
were analyzed, and a baseline was defined, the program plan and projected 
costs were generated, and a risk assessment was made. To perform these 
design studies it was necessary to postulate or extrapolate extended- 
performance characteristics in several technology areas. Consequently, 
studies were conducted in four technology areas to confirm the validity 
of the design approaches that were considered. 

The thrust systems that were considered are based on the electric 
propulsion technology that NASA LeRC has been developing for over a 
decade. The technical baseline for this application is the most recent 
operational engineering model thruster (EMT), the 900-series 30-cm mercury 
ion EMT. This thruster is a scaled-up version of the 15-cm thruster 
developed and flight tested during the 1960-1969 period for the SERT II 
program. The EMT operates at a 3-kW power level with a specific impulse 
of 3,000 sec. By making minor modifications in the existing thruster 
design, extended performance at approximately 6-kW power level, 4,800-sec 
specific impulse, and 15,000-hr pre-wearout life (as required for a 
Halley's coitjet mission) was believed to be achievable at a low technical 
risk.- This supposition was evaluated as part of this study. 
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In addition to the extended performance thruster, the key elements 
of the thrust system for this extended-performance application are the 
power management and control (PMaC) subsystem, the gimbal system, pro- 
pellant storage and distribution system, thermal control, and supporting 
structure. The background of extensive development in power-processing 
technology for mercury ion thrusters and technology developments in the 
other areas were the basis for the high level of confidence that the 
required extended performance levels could be achieved. 

B. SCOPE 

The scope of this study included’: the development of conceptual 

designs for various candidate systems; the selection, definition, and 
evaluation of a baseline design concept and its critical interfaces; an 
evaluation of the sensitivity of the baseline design to critical data 
base and design parameters; the generation of a development program plan 
for the baseline concept; estimation of costs and fiscal year funding 
requirements; fabrication of a demonstration scale model; and the conduct 
of supporting technology studies (including fabrication and testing of 
critical hardware components) to estimate the physical and electrical 
performance and to provide a baseline for subsequent work. 

The design characteristics, program plan, and costs of the baseline 
system, were defined in parallel with the supporting technology effort. 
Design definition was carried out in two consecutive phases: 

• Phase 1: Definition and comparison of alternate 

configurations, leading to baseline selection 

• Phase 2: Design definition and evaluation of the 

baseline configuration, culminating in the generation 
of a program plan and cost estimates. 

The evaluation of supporting thruster technology consisted of three 
essentially independent investigations to confirm the critical assump- 
tions made during the design study in the following areas: 

• Operation and documentation of a modified 900-series 
EMT at the beam voltage and power levels required for 
a Halley's comet mission 
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• Design and evaluation of thruster propellant electrical 
isolators for operation at beam voltages up to 4,000 V 

• Development of a model for the charge-exchange plasma 
efflux from the ion beam exhaust plume (for analysis 
of interactions between the thrust system and a solar 
array operated at high voltage). 

The program schedule, which required that the results of the completed 
study be reported during July 1977, determined the approach and scope of 
the thrust system design study and of the supporting studies. Much of the 
design study was therefore limited to conceptual definition of the key 
characteristics of subsystems and components, with just enough detail 
included to permit realistic estimates of mass, power dissipation, and 
reliability. Similarly, we sought to confirm the technology approaches 
in question by using straightforward engineering extrapolations of 
existing data or hardware designs, and the scope of the study was quite 
limited. 

The fewest possible modifications were made to the 900-series EMT 
and evaluated. The high-voltage isolator design study was based on 
already proven approaches; and the charge-exchange plasma model was 
based on data available from previous ion beam measurement and theoretical 
analyses. Since the intent was to evaluate feasibility, conclusive 
verification of any specific design is not implied.. Although the 
thruster technology areas do relate directly to the design studies 
described in Volumes I, II, and III, the results are presented here, in 
a separate volume, because those areas were investigated independently. 
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SECTION 2 


MODIFICATION AND EVALUATION OF THE 900-SERIES 30-cm EMT 
FOR EXTENDED-PERFORMANCE OPERATION 

Recent developments 1n the technology of thin solar cells, innovative 
solar array designs, and advanced technology ion thruster development 
have made feasible energetic missions such as the Halley's comet rendez- 
vous. These high-energy-requirement missions, which require propulsion 
system specific masses of 10 to 25 kg/kW, would use the launch capabili- 
ties of the Space Shuttle's intermediate upper stage (lUS) combined with 
solar electric propulsion employing 30-cm diameter electron-bombardment 
mercury ion thrusters. The 900-series 30-cm EMT, developed for primary 
propulsion applications, is now close to achieving flight hardware status. 

The EMT has been successfully operated and performance documented over a 
throttling range of 4:1;^ it has a demonstrated lifetime of 10,000 hr (Ref. 2). 
A mathematical model of this thruster was recently developed that describes 
the performance parameters and thrust correction factors with respect to 
the independent variables of beam current and beam voltage. 

A. DEFINITION OF THE EXTENDED PERFORMANCE OPERATING RANGE AND 

THRUSTER REQUIREMENTS 

The design conditions of the 30-cm EMT were selected early in its 
development program to satisfy the requirements of a typical baseline . 
set of planetary and Earth orbital missions using the solar array and 
power conditioning technology (resulting in a propulsion system specific 
mass, a, of 50 kg/kW) that was available at that time. The specific 
impulse design point of the 30-cm EMT is 3,000 sec, which corresponds to 
a beam voltage of 1,100 V. ' At this voltage, a thrust level of 128 mN is 
produced when operating at the beam current design point of 2 A. 

The propulsion system specific mass of 10 to 25 kg/kW, a solar array 
power of 50 to 100 kW, and other Halley's comet mission requirements 
dictate the use of a thruster that has a higher power and specific 
impulse than does the present 900-series 30-cm EMT. Figure 1 presents 
thruster operating envelopes for the various pov/er systems considered in 
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Figure 1. Thruster operating envelopes for various PMaC approaches 
proposed for the Halley's comet rendezvous. 
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the design study for the Halley's comet mission. Also shown is the 
operating range over which the 900-series EMT has been performance docu- 
mented and mathematically modeled. The 900-series EfIT does have the 
capability to operate below the beam voltage range shown in Figure 1. 
Comparing these envelopes illustrates how far some of the PMaC approaches 
considered depart from the basic EMT operating range. To operate the 
EMT in this extended performance range would require at least two 
modifications to the current design: 

e Propellant isolator redesign or modification to permit 
thruster operation at high levels of beam voltage. 

• An increase in the accelerator system gap to reflect the 
lower perveance requirements of the extended-performance 
thruster and to provide increased voltage standoff 
capability. 

In addition, several factors peculiar to high-voltage thruster opera- 
tion must be considered. These include the adequacy of wire insulation, 
protection of insulator surfaces, and close examination of threaded 
parts and sharp edges on insulator shields. 

Preliminary verification and documentation of thruster operation 
at high power levels was accomplished under NASA contract NAS 3-19703 
(Ref. 3). During that program, high beam current performance tests 
were conducted over a beam current range of 2 to 4 A with a nominal 
beam voltage of 1,100 V. Performance at high beam voltage levels was 
documented over a voltage range from 2 to 5 kV while maintaining a 2-A 
beam current. The test performance results were in agreement with pre- 
dictions. The major conclusions from these preliminary studies were that 

• High-voltage operation has no discernible effect on 
discharge chamber wear. rates, but requires adequate 
electrical insulation. 

• ' High-current operation increases the discharge chamber 

wear rate (and shortens thruster lifetime as a 
consequence) . 

Although those preliminary tests verified that the 30-cm EMT could be 
operated in the extended-performance range with relatively few modifica- 
tions, performance data was still needed for a much broader operating 
range, as indicated in Figure 1. This documentation was obtained 
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under the' test program described here. The primary objectives of the 
investigation were to; 

• Verify that the modified EMT can be operated throughout 
the performance envelopes representative of the Halley's 
comet rendezvous mission 

• Assess the applicability of the EMT performance model in 
the increased beam current and beam voltage ranges. 

Experimental evaluation of the extended- performance thruster included:’ 

• Performance documentation 

• Measurement of propellant utilization and thrust correction 
factors 

• Measurement of beam cone angle 

o Determination of discharge chamber and screen grid erosion 
rates 

• Measurement of various thruster component temperatures. 

B. THRUSTER AND TEST FACILITY DESCRIPTIONS 

To operate the 30- cm EMT at beam voltage levels representative 
of the Halley's comet rendezvous mission required that the thruster 
modifications discussed in the previous section be made. The modified 
thruster used during the test program, the vacuum test facility, and 
laboratory power system are described in Sections 2.B.1 and 2.B.2. 

1 . Thruster 

Figure 2 is a photograph of the extended-performance thruster 
(SN 801) used during the test program. This thruster is an 800-series 
thruster modified to include all the components from the 900 series 
that were expected to affect thruster lifetime and performance charac- 
teristics. These modifications included: 

• Flake conthol treatment of various discharge chamber 
surfaces 
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Figure 2. Photograph of thruster SN 801 modified for 
extended performance testing. 
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• Erosion control treatment of selected discharge-chamber 
and neutralizer surfaces 

• 400-series thruster design baffle mount 

t Swaged-type cathode heaters 

• Impregnated porous tungsten cathode inserts 

• A separable-type (cathode from vaporizer) neutralizer 
assembly. 

The flake control surface treatment consisted of a wire mesh 
covering on the anode and the upstream discharge chamber boundary. 

Erosion control treatment involved installing tantalum cladding on the 
cathode pole piece and baffle and applying thin graphite strips to the 
neutralizer housing. In addition to these modifications, insulated 
restraints were installed at various locations to prevent wires from 
contacting thruster structural members. 

As discussed earlier, the reduced perveance requirements of the 
extended-performance thruster allows the grid separation to be 
increased, thereby providing an improved voltage standoff capability. 

The EMT perveance and the perveance requirements of the extended- 
performance thruster are indicated in Figure 3. The perveance require- 
ments of both thrusters are dictated by the minimum total extraction volt- 
age available at the 2-A beam current level; Figure 3 shows the minimum 
voltage to be about 700 V higher for the extended-performance thruster 
configurations. The increase in grid spacing corresponding to the line 
labelled "required perveance" in Figure 3 can be calculated by using the 
EMT grid separation of 0.5 mm (0.020 in.) and the total extraction volt- 
ages at the 2-A beam current level, and by assuming that the current is 
governed by space-charge-limited flow between planar electrodes. The 
grid separation calculated in this manner is 0.76 mm (0.030 in.). This is 
the nominal spacing used during the test program and represents a 25% 
reduction in the grid separation used in the preliminary high-voltage 

3 

thruster tests. Figure 4 presents the variation in measured grid separa- 
tion for ion optics assembly SN 817, which was used during the present 
program. Under atmospheric conditions, applying grid voltages greater 
than 2,500 V produced arcing in the minimum grid separation region 
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Figure 3. Thruster perveance requirements for different 
PMaC approaches. 
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ACCELERATOR ELECTRODE 


APERTURE DIAMETER 1.52 

SPACING CTR/CTR 2.21 

THICKNESS 0.51 

SCREEN ELECTRODE 

APERTURE DIAMETER 1.91 

SPACING CTR/CTR 2.21 

THICKNESS 0.38 

REDUCTION 0.4 % 



Figure 4. Electrode spacing chart for EMT grid design with nominal 

interelectrode spacing increased from 0.5 to 0.76 mm (dimen- 
sions are in mm except as noted). 
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indicated in the figure. However, before the thruster test program was 
begun, voltage standoff capability under high-vacuum conditions was 
demonstrated with applied voltages as high as 4,400 V. 

One of the first tests conducted with the extended-performance 
thruster .was to verify that the perveance and voltage standoff capabil- 
ity of the re-spaced ion optics were adequate for the extended perform- 
ance operating envelopes shown in Figure 3. The figure shov/s that the 
perveance requirement is set by the operating point corresponding to 
beam conditions of 2 A and 1600 V (an accelerator voltage of 500 V was 
used in all tests). Perveance data was obtained at these conditions 
by reducing beam voltage while keeping discharge power and accelerator 
voltage constant. Figure 5 presents the measured variation of accel- 
erator and beam currents with respect to total extraction voltage. The 
perveance line constructed from the measured variation in accelerator 
current, the EMT perveance line, and the required extended-performance 
thruster perveance lines are also shown in Figure 5. This data indicates 
that the perveance of the EMT optics, re-spaced to a separation of 0.76 mm, 
is more than adequate for operation in the performance envelopes charac- 
teristic of the Halley's comet mission. Subsequent testing of the 
extended-performance thruster verified the adequacy of high-voltage 
insulation and ion optics voltage standoff capability up to the 10-kW 
power level that corresponds to beam conditions of 2.5 A and 3,900 V. 

2., Test Facility 

The Hughes Research Laboratories (HRL) 2.7 m (9 ft) diameter 
vacuum facility was used for all thruster testing. The tank is cryo- 
pumped and has maintained ambient chamber pressure on the order of 

r 

10~ Torr for ion beam power levels up to 10 kW and beam currents up 
to 4 A. Power to the various thruster components was provided by 
general-purpose 1 aboratory-type supplies. These power supplies were 
adequate for thruster operation throughout the extended performance 
envelopes indicated in Figure 1. A block diagram of the thruster power 
system is presented in Figure 6. 
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Figure 5. Perveance data for EMT grid design with interelectrode spacing 
increased to 0.76 mm. 
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Figure 6. Block diagram of the power system of the extended- performance thruster 











The design and testing of high-voltage propellant isolators was 
done concurrently with the thruster testing program. This required 
short-circuiting the EMT cathode and main isolators to allow these 
propellant feed systems to float at thruster body potential. Conse- 
quently, the load cell flow measurement technique normally used for 
30-cm thruster testing and documentation was not available for use 
under this program because the main and cathode feed systems were 
floated as high as 3.9 kV, which is beyond the isolation capability of 
the load cell instrumentation. Instead, the flow rates were determined 
by timing the fall of liquid mercury columns in calibrated glass 
pipettes. 

In addition to providing power for the various thruster components, 
the power system assembled for this test program allowed certain 
operating parameters (such as beam current) to be held at a constant 
value. This was done by using current- or voltage-limited power 
supplies and closed-loop vaporizer control systems. Beam current was 
maintained at the desired value by controlling the main vaporizer 
current; the discharge and neutralizer keeper voltages were maintained 
at selected values by closed-loop control of the cathode and neutralizer 
vaporizer currents. The cathode and main vaporizer control loops were 
considered necessary to ensure that the data acquired would be accurate 
and repeatable for relatively short test times (< 2 hr) and to main- 
tain test parameters constant over longer periods (up to 30 hr) for 
tests using the f x t momentum analyzer probe and for erosion rate 
measurements. 

C. SELECTION OF THRUSTER OPERATING PARAMETERS 

The values of beam current and voltage selected for thruster 
performance documentation tests are representative of the operating 
envelopes required for the various PMaC approaches proposed for the 
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Halley's comet rendezvous mission. These envelopes and the operating 
points selected for the performance documentation tests are presented 
in Figure 7. The discharge parameters corresponding to the selected 
operating points are presented in Table 1 . A discharge power loss of 
200 eV/ion for these operating points was chosen because it would 
yield good propellant utilization at acceptable power levels. With 
the exception of some limited parameter variation testing, the discharge 
voltage was held constant at 36 V, which is a reasonable compromise 
between thruster efficiency and lifetime considerations. 

An accelerator grid voltage of 500 V was chosen to eliminate 
electron backstreaming without causing excessive beam divergence. 

This choice was made on the basis of preliminary high-voltage thruster 
tests conducted under NASA Contract NAS 3-19703. These tests had 
revealed that thruster operation at beam voltages higher than 2 kV 
causes electron backstreaming if the absolute value of accelerator volt- 
age is less than 400 V. Subsequent testing under the current program of 
the extended-performance thruster confirmed that the minimum accelerator 
grid voltage necessary to prevent electron backstreaming was 400 V. 
Cathode and neutralizer keeper currents of 0.8 and 2.0 A, respectively, 
were selected to be consistent with the current EMT operating conditions. 
The magnetic baffle current was set at either 2.2 or 2.7 A (approxi- 
mately). The lower value corresponds to the minimum cathode flow rate 
required to prevent inefficient "low-mode" discharge-chamber operation; 
the higher value was selected to allow evaluating the effects of 
cathode flow rate on the thruster performance parameters. 

D. TEST PROCEDURE 

The thruster was operated at the desiVed conditions- for approxi- 
mately 3 hr after startup. By this time the starting transients had 
decayed, the tank pressure had reached the low 10 ® Torr range, and 
the accelerator grid current was typically less than 0.2% of the beam 
current. Propellant flow rate measurements were then started and 
continued for approximately 2 hr. These measurements consisted of 
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Figure 7. Thruster operating points selected for performance 
documentation. 
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Table 1. Operating Parameter Range for the 
Extended Performance Thruster 


Beam 

Discharge 


''b- 


''d’ 

A 

kV 

A 

V 

2.5 

3.9 

13.9 

36 

2.5 

3.0 

13.9 

36 

2.5 

2.5 

13.9 

36 

2.0 

3.0 

11.4 

35 

2.0 

3.0 

11.1 

36 

2.0 

2.5 

11,1 

36 

2.0 

1 .6 

n.i 

36 

1.5 

3.0 

8.3 

36 

1.5 

2.5 

8.3 

36 

1.5 

1.6 

8.3 

36 

1 .0 

3.9 

5.6 

36 

1.0 

3.0 

5.6 

36 

1.0 

2.5 

5.6 

36 

1.0 

1.6 

5.6 

36 

a 

200 eV/ion discharge loss 
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recording the height of a liquid mercury column at 10 min time intervals. 
The slope of the volume-versus-time curve plotted from this data was 
used to determine propellant flow rate. The plots were used to provide 
a visual indication of linearity and, hence, of steady-state flow 
conditions. After sampling flow data for approximately 2 hr, the 
flow rate for each vaporizer was calculated using linear least-squares 
regression. The correlation coefficient, which indicates the constancy 
of the flow rate during the data sampling period, was typically greater 
than 0.999. 

Electrical measurements and tank pressure readings were also 
recorded at 10-mi n time intervals. Currents and voltages critical 
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to thruster performance calculations (discharge, beam, and' accelerator) 
were measured using ±1% accuracy mirrored-scale panel meters. All other 
electrical quantities were measured using ±3% accuracy meters. 

Beam diagnostic' data was obtained using an f x ^ momentum analyzer 
probe. This instrument was used to measure the positively charged 
constituents of the ion beam and to measure thrust losses due to beamlet 
divergence and the presence of multiply charged ions. Beam profiles were 
recorded using a single Faraday cup probe that was swept through the ion 
beam at various axial locations. 

E. EXPERIMENTAL RESULTS 

The thruster test program consisted of two basic tasks: performance 

documentation and assessment of thruster lifetime. Performance docu- 
mentation involved determining thruster performance parameters (such as 
electrical and propellant utilization efficiencies) and measuring the 
efficiency correction factors (which correct for losses caused by beam 
divergence and the presence of multiply charged ions). The thruster life- 
time investigation, which consisted of three types of tests, was con- 
ducted to enable assessing the effects of high-power thruster operation 
on thrust system lifetime. These tests included measurements of discharge 
chamber and screen grid erosion rates, beam cone angle,* and critical 
thruster component temperatures.' 

1 . Performance Documentation 

Performance testing consisted of obtaining flow rate, electrical, 
and beam diagnostic measurements at each of the beam conditions indicated 
in Figure 7. The performance parameters were calculated using the elec- 
trical measurements indicated in Figure 6 and the following expressions 
for total power 

V<''b- V^<-'d+’b + ‘a)''d (1) 

^VA^^^raisc • 

*Cone angle is the angle of a conically shaped region that contains 
95% of the primary ion beam ion current. 
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electrical efficiency 


= 


'b''b 


e P, ■ 


propellant utilization efficiency 


n.. = 


u m, 


B 

TOTAL 


and thruster efficiency 




( 2 ) 


(3) 


(4) 


where total propellant flow rate expressed in equivalent 

amperes, and represents the miscellaneous thruster input power 

(which consists of magnetic baffle, keeper, and vaporizer heater powers). 
The thrust correction required to account for the presence of doubly 
charged ions, a, and the thrust loss correction due to beamlet divergence, 
F^, were calculated from the f x ^ probe measurements in the manner 
described in Ref. 4. 

The calculated performance data is presented in Table 2. This per- 
formance data is compared with the values calculated using the EMT 
mathematical model in the following section. There are trends evident 
in the performance data that are not accounted for in the model, 
however, and these effects are discussed below. 


a. Beam Voltage Effects 

The propellant utilization data (Table 2) generally indicates 
that utilization increases with increasing beam voltage. This trend, 

3 

which was also observed in the preliminary high-power thruster tests, 
probably is a result of plasma sheath movement within the screen grid 
apertures. As beam voltage is increased with beam current held 
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Table 2. Thruster Performance Data for Extended-Performance 

Operations 


Controlled Parameters^ 



Measured 

Parameters 



Beam Diagnostics 

B 

Vb> 

V 

A 

Vd. 

V 

Imb* 

A 

lA. 

toA 

A 

^c* 

A 

®T0TAL 

A 

a,b 
^ TOTAL, 

b 

b 

Hu- 

X 

b 

nt5 

z 

i+. 

% 

a 

e 


“ 


1.0 

1.6 

5.6 

36 

2.3 

m 

1.147 

0.082 

1.259 

1,877 

85.2 

79.4 

62.2 

8.4 

0.981 

0.961 

0.977 

0.919 

0.14 

1.0 

1.6 

5.6 

36 

3.0 

1.8 

1.096 

0.106 

1.232 

1,886 

87.8 

81.2^ 

63.0 

8.9 

0.980 

0.959 

0.976 

0,915 

0-21 

1,0 

2.5 

5.6 

36 

2.3 


1.156 

0.086 

1.272 

2,778 

90.0 

78.6 

65.0 

7.9 

0.982 

0.963 

0.976 

0.919 

0.15 

1.0 

2.5 

5.6 

36 

3.0 

H 

1.109 

0.106 

1.245 

2,787 

89.7 

80.3 

65.6 

8.2 

0,981 

0.962 

0.973 

0.911 

0.19 

1.0 

3.0 

5.6 

36 

2.3 

H 

1.140 

0.088 

1.250 

3,279 

91.5 

80.0 

66.8 

8.0 

0.982 

0.963 

0.973 

0.913 

0.16 

1.0 

3.6 

5.6 

36 

2.8 

H 

1.075 

0.115 

1.220 

3,285 

91.3 

82.0 

68.2 

9.2 

0.979 

0.958 

0.975 

0.911 

0.15 

1.0 

3.9*^ 

5.6 

36 

2.6 

1.8 

1.032 

0.153 

1.215 

4,182 

93.3 

82.3 

70.0 

9.2 

0.979 

0.958 

0.975 

0.911 

0.15 

1.3 

1.6 

8.3 

36 

2.2 

3.1 

1.651 

0.074 

1.755 

2,778 

86.4 

85.5 

67.0 

11.9 

0.974 

0.997 

0.976 

0.907 

0.27 

1.5 

1.6 

8.3 

36 

2.3 

3.4 

1.630 

0.075 

1.735 

2,779 

86.4 

86.5 

67.9 

11.6 

0.974 

0.973 

0.979 

0.969 

0.26 

1.5 

1,6 

8.3 

36 

3.0 

2.7 

1.566 

0.094 

1.690 

2,771 

86.6 

88.8 

69.2 

13.2 

0.971 

0.942 

0.977 

0-900 

0.35 

1.5 

2.5 

8.3 

36 

2.2 

2.7 

1.632 

0.079 

1.741 

4,129 

90.8 

86.2 

70.9 

12.0 

0.973 

0.946 

0.978 

0.906 

0.27 

1.5 

2.5 

8.3 

36 

3.0 

2.6 

1.531 

0.113 

1.679 

4,122 

91.0 

89.3 

73.9 

12.1 

0.973 

0.946 

0.980 

0.909 

0.30 

1.5 

3.0 

8.3 

35 

2.2 

2.5 

1.611 

0.080 

1.721 

4,879 

92.2 

87.2 

72.8 

11 4 

0 974 

0.949 

0.977 

0.906 

0.30 

1.5 

3.0 

8.3 

36 

2.7 

1.9 

1.585 

0.085 

1.700 

4,884 

92.1 

88.2 

72.9 

13.4 

0.970 

0.941 

0.977 

0.898 

0.31 

1.5 

3.0 

8.3 

36 

2,8 

1.9 

1.565 

0.101 

1.696 

4.225 

92.1 

88.4 

72.7 

11.5 

0.974 

0.949 

0.977 

0.893 

0.28 

1.5 

3-0 

8.3 

36 

3.0 

2.7 

1.501 

0.124 

1.655 

4,873 

92.3 

90.6 

75.3 

11.8 

0.974 

0.977 

0.974 

0.900 

0.31 

2.0 

1.6 

11.1 

36 

2.3 

4.3 

2.100 

0.063 

2.193 

3,683 

86.9 

91.2 

70.3 

18.3 

0.961 

0.923 

0.980 

0.887 

0.43 

2.0 

1.6 

11.1 

36 

2.3 

4.5 

2.088 

0.058 

2.176 

3,684 

86.9 

91.9 

71.2 

16.4 

0.965 

0.930 

0.978 

0.891 

0.72 

2.0 

1.6 

11.1 

36 

2.8' 

3.5 

2.039 

0.093 

2.162 

3,691 

86.7' 

92.5 

71.1 

18.4 

0.961 

0.923 

0.980. 

0.887 

V 0.49 

2.0 

2.5 

11.1 

36 

2.3 

3.6 

2.094 

0.069 

2.193 

5,483 

91.2 

91.2 

79.1 

16.6 

0.965 

0.929 

0.978 

0.891 

0.51 

2.0 

3.0 

11.4 

35 

2.2 

3.4 

2.079 

0.074 

2.183 

6,484 

92.5 

91.6 

75.7 

15.4 

0.967 

0.933 

0.978 

0,894- 

0.77 

2.0 

3.0 

11.1 

36 

2,3 

3.3 

2,048 

0.073 

2.151 

6,484 

92.5 

93.0 

76.1 

17.3 

0,963 

0.926 

0.977 

0.885 

0.50 

2.0 

'3.0- 

11.1 

36 

2.7 

2.6 

2.016 

0 lOS 

2.154 

6,491 

92.4 

92.9 

75.5 

19.5 

0.959 

0.918 

0.978 

0.880 

0.51 

2.0 

3.0 

11.1 

36 

2.7 

2.8 

2.001* 

0.112 

2.143 

6,492 

92.4 

93.3 

76.4 

17.7 

0.962 

0.925 

0.489 

0.885 

0.48 

2.0 

3.0 

11.1 

36 

2.7 

2.4 

1.977 

0.113 

2.120 

6,490 

92.4 

94.3 

76.8 

19.0 

0.960 

0.920 

0.978 

0.881 

0.50 

2.5 

2.5 

13.9 

35 

2.2 


2.467 

0.078 

2.575 

6,842 

91.4 

97.1 

76.6 

24.5 

0.951 

0,902 

0.977 

0.863 

0.74 

2.5 

3.0 

13.9 

36 

2.2 

B 

2,451 

0.088 

2,569 

8,092 

92.7 

97.3 

77.4 

24.4 

0.951 

0.902 

0.974 

0.858 

0.70 

2.5 

3.9'^ 

13.9 

36 

2.7 

B 

2.336 

0.119 

2.485 

10,347 

94.2 

100.6 

81-3 

24.4 

0.951 

0.912 

0.974 

0.858 

0.70 


^Accelerator voltage was 500 V. 

*^Assumes = 13.8 V, = 1.8 A, ~ 10 V, neutralizer flow rate = 0.030 A. 
^Beam diagnostic data was obtained at 3 kV beam voltage. 
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constant, the plasma sheath moves upstream, which reduces ion loss to 
the screen aperture walls and increases the effective open area of the 
screen grid. These effects permit discharge operation at. a lower 
plasma density, which implies lower neutral density and increased 
propellant utilization at a fixed beam current level. 

b. Magnetic Baffle Current Effects 

Reducing the magnetic baffle current and the cathode flow rate 
proportionately generally improved the performance of the discharge 
chamber. Improved performance at low cathode flow conditions is 

apparently a consequence of fewer neutrals passing through the baffle 
aperture region. This permits a larger fraction of the electrons drawn 
from the cathode discharge plasma to attain an energy corresponding to 
the potential difference between the cathode and main plasmas, which 
results in an increase in the ionization cross section and more effi- 
cient thruster operation. However, a lower limit on cathode flow rate 
exists since very low cathode flows generally result in an inefficient 
mode of discharge operation (sometimes referred to as "low mode"), 
characterized by poor propellant utilization efficiency. 

Table 2 presents thruster performance data for both low and high 
magnetic baffle currents (approximately 2.2 and 2.7 A), corresponding 
to cathode flow rates in the range from 60 to 150 mA. This data shows 
a consistent trend towards more efficient thruster operation at higher 
cathode flow rates. This is shown by the increase in propellant 
utilization and total thruster efficiency which accompanied increased 
cathode flows. The increase in propellant utilization at constant 
beam current generally reduced accelerator drain current and (presum- 
ably) erosion of the accelerator grid caused by charge-exchange ion 
impingement. 

The reason for the performance improvement that was observed at 
higher cathode flow rates is not well understood. A complete under- 
standing would require acquiring plasma diagnostic and magnetic field 
geometry information, and this was beyond the scope of this program. 
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Lacking this data, all that can be said is that the extended- performance 
thruster may- be operationally different from the unmodified 900rser.ies 
EMT in this respect. If such a difference exists, it would need to 
be determined whether it is caused by plasma, geometric, or magnetic- 
field effects. Efficient and stable thruster operation at higher cathode 
flow rates is considered a desirable characteristic since it .would 
provide a wider margin for flow variation above the low-mode lower limit 
on cathode flow rate. 

c. High Power Operation 

The performance data in Table 2 indicates that electrical 

efficiency is higher when operating at higher power levels. This 

occurs because the discharge and miscellaneous fixed power losses 

become smaller fractions of the total thruster power as the beam 

power is increased. As discussed earlier, operation at increased 

beam voltages with beam current held constant increases propellant 

utilization efficiency. Also, at a fixed beam voltage, the propellant 

utilization efficiency increases with beam current level. The increase 

in propellant utilization with beam current is expected since thruster 

operation at the 200 eV/ion discharge power level generally results 

5 

in a constant neutral loss rate. 

d. Discharge Voltage Reduction 

Lowering the discharge voltage level to 35 V reduces 
propellant utilization and the ratio of double-to-s ingle ions. The 
combination of reduced plasma potential (ion energy) and double ion 
density achieved through discharge voltage reductions is considered 
an effective means for increasing thruster lifetime. However, the lowered 
discharge voltage reduces the cross sections for the formation of both 
single and double ions. This results in a loss in propellant utilization 
efficiency. Acceptable propellant utilization efficiency may be main- 
tained only if the neutral density is increased to offset any reduction 
ion ion production rate. 
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An acceptable neutral loss rate can be maintained at the increased 
neutral density by reducing the open area of the accelerator grid. The 
use of small-hole accelerator grid (SHAG) optics for obtaining good 
thruster performance at low discharge voltage levels was demonstrated 

O 

under a previous program. No effort was made during the current 
program to verify the effectiveness of SHAG optics in improving ’thruster 
lifetime. 

2. Thruster Lifetime 

Three types of tests were conducted to assess the lifetime capabil- 
ities of the extended-performance thruster: 

• Beam divergence angle measurements for use in assessing 
the compatibility of the thruster and the solar array 
and defining this interface requirement. 

• Temperature measurements of various thruster components 
for determining the effects of high-power thruster 
operation on local temperatures near critical areas 
(such as propellant isolators, wire insulation, and braze 
joints) . 

• Erosion rate measurements on critical thruster components 
(such as the screen and accelerator grids, pole piece, 
baffle, and baffle supports). 

The results of these tests are discussed separately below. 

a. Beam Divergence 

The primary ion beam cone angle was determined by sweeping 
a Faraday probe through the beam while simultaneously recording probe 
position and collector current. These measurements were made at 
axial locations 10 and 20 cm downstream from the accelerator grid; 
the beam current that was obtained by integrating the current density 
profiles measured with the Faraday probe generally agreed with the 
metered beam current to within ±5%. 

The radial position at which the current density drops to zero 
can be used to define the beam edge, although this method of defining 
the beam edge is considered somewhat inaccurate because of the 
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detection threshold and finite collection area of the- Faraday probe. 

An alternate and perhaps more meaningful procedure is. to define the 
beam edge as that radius which contains a large fraction (say 955^) of 
the total primary ion beam. This approach reduces any error introduced 
by the detection threshold of the probe, since the "tail" of the current 
density profile contributes a relatively small amount to the beam 
current integral . 

Beam cone angles obtained from the Faraday probe measurements are 
presented in Figure 8 as a function of beam current. The angles were 
calculated using the beam radius at a known axial location, and a beam 
radius at the thruster exit plane of 14.5 cm. The variation of beam 
cone angle with beam voltage was found to be negligible, which is not 
surprising since the net-to-total beam voltage ratio, R, remained 
relatively constant. Figure 8 indicates that a 20° cone contains 95% 
of the primary beam ions and that this value is essentially independent 
of beam conditions. The beam cone angles calculated from the radius 
at which the Faraday probe signal drops below measurement threshold (less 
than 6 yA/cm ) is also shown in Figure 8; this data suggests a maximum 
beam cone angle of approximately 40°. The absolute beam edge determined 
in this manner is considered optimistic, since recording the Faraday 
probe output on a more sensitive current indicator would tend to increase 
the calculated radius at which the signal drops below the measurement 
threshold. Based on probe measurements made using more sensitive current 
scales, the error in absolute beam edge is estimated to be within -i-8°. 

An alternate and more accurate procedure for determining the 
absolute beam edge was accomplished through the use of alternating-film 
erosion monitors.* These monitors were attached to the neutralizer 
housing, from which an unobstructed view of the primary ion beam edge 
was provided. This technique accurately defines the beam edge by 


'J'C 

A complete description of the alternating-film erosion monitor 
technique is presented later in this section. 
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Figure 8. Measured beam cone angles. 


27 




determining the axial location on the neutralizer housing upstream 
of which no sputtering damage occurs. Beam cone angles determined in 
this manner are presented in Table 3. This data indicates that the 
beam cone-angle is k 50° and is independent of beam current. An 
increase in the beam cone angle of approximately 1° was observed at 
the lower beam voltage, and this increase in beam divergence was 
apparently due to the reduction in the net-to-total beam voltage 
ratio, R. 

Beam divergence angles determined using the erosion monitors have 
an estimated accuracy of ±1°, with the bulk of the error caused by 
errors in geometric measurements. Comparison of the absolute beam 
cone angles determined by the Faraday probe and erosion monitor tech- 
niques indicates agreement to within 2° when the error in the Faraday 
probe method is accounted for. 

b. Temperature Measurements 

Thruster component temperatures were measured using Chromel- 
Alumel thermocouples attached at the critical locations indicated in 
Figure 9. After several hours of thruster operation at the 6.5-kW 
power level (2-A, 3000-V beam conditions), the indicated temperatures 
were recorded using a digital thermocouple meter. These measurements 

are in good agreement w.ith thermal data obtained during preliminary 

. 3 

high-power thruster testing at the 5-kW beam power level. The 

temperatures recorded at the 6.5-kW power level are only 40 to 50°C 

higher than those corresponding to EMT operation at 2.5 kW and are 

considered to be within the acceptable temperature range of these 

components . 

c. Erosion Rate Measurements 

The erosion rates of critical thruster components were 
measured at various thruster operating conditions using the multilayer 

3 

erosion monitor technique developed under previous NASA and COMSAT 
contracts. The thruster components selected for erosion rate 
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Figure 9. Schematic of thruster discharge chamber showing thermocouple 
locations and measured temperatures corresponding to 2-A, 
3-kV beam conditions. 
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Table 3. Beam Cone Angles Determined from Alternating- Film 
Erosion Monitor Measurements 


A 

''b> 

kV 

9, 

deg 

R® 

1.0 

3.0 

50.5 

0.86 

1.5 

3.0 

50.4 

0.86 ■ 

2.0 

3.0 

50.1 

0.86 

2.0 

1.6 

51.4 

0.76 

E ratio of the net-to-total accelerating voltages 
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measurements included the screen and accelerator grids, baffle and 
support, cathode pole piece, and neutralizer housing. 

d. Erosion Monitor Description 

The multilayer erosion monitors, which are prepared in the 
Sputtering Processes Laboratory at HRL, consist of alternating layers 
of color-contrasting materials (such as tantalum and copper, or molybdenum 
and copper) that are sputter-deposited on a tantalum or stainless-steel 
substrate. Layer thickness is carefully controlled during the applica- 
tion process by maintaining the sputtering ion beam voltage and current 
constant and by precisely timing the deposition. A small piece of 
polished material sputter-deposited during preparation of the multi- 
layer monitor material is used for calibrating layer thickness using 
a surface profilometer. Typical layer thickness for erosion monitors 
produced and calibrated in this way is 60 ±0.6 nm for the usual deposi- 
tion control conditions. 

Multilayer monitors used inside the discharge chamber had 8 tan- 
talum or molybdenum layers and 7 copper layers; the first layer (either 
tantalum or molybdenum) applied to the substrate was double thickness. 
Because copper erodes about 20 to 40 times faster than either molyb- 
denum or tantalum, the copper layers are factored into the analysis 
as only a portion of a layer. Therefore, the total equivalent thickness 


30 













of a 15 layer (8 Ta or Mo plus 7 Cu layers) monitor having 60 nm layer 
thickness is approximately 566 nm of molybdenum or 553 nm of tantalum. 

e. Test Procedure 

Small pieces of the thin- film erosion monitor material were 
spot-welded to the thruster components selected for erosion measure- 
ments. A portion of the monitors mounted on all thruster components 
except the screen and accelerator grids was masked to aid in determin- 
ing the number of layers that had been eroded during thruster operation. 
Those monitors mounted on the screen and accelerator grids were per- 
forated, with the monitor hole pattern matching that of the screen 
grid. Photographs to show typical examples of the erosion monitors 
mounted on various thruster components are presented in Figure 10. A 
summary of the location, type, and thickness of the erosion monitors 
used during these tests is given in Table 4. 

Several precautions were taken during the thruster erosion monitor 
testing to ensure valid data acquisition. Among these were careful 
control of the discharge chamber operating conditions by means of 
vaporizer control loops, and checks on the double ion content of the 
discharge chamber using the E x B probe. Test chamber vacuum condi- 
tions are crucial in obtaining valid erosion-rate measurements since 
the presence of adsorbed gases and surface contamination (from poorly 
controlled vacuum conditions) could substantially reduce the surface 
sputtering yield. Several procedures were followed during these tests 
to ensure accurate erosion rate data by carefully monitoring and con- 
trolling the vacuum chamber pressure. The control procedures 
included: 

• Venting the vacuum tank with gaseous nitrogen before 
removing or installing the thruster 

• Sufficient pumping after thruster installation to achieve 
chamber pressures in the low 10“^ Torr range without 
cryo-pumping. 
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Figure 10. Photographs of multilayer erosions monitors mounted on 
internal thruster surfaces. (From Ref. 3.) 
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Table 4. Erosion Monitor Locations and Layer Description 


Thruster 

Location 

Monitor 
Materi al 

Layer 
Thickness , 
nm 

Number of 
Layers 

Screen and 
accelerator grid 

Mo and Cu 

60 

8 Mo + 7 Cu 

Baffle, baffle 
support, cathode 
pole piece 

Ta and Cu 

60 

8 Ta + 7 Cu 

Neutral izer 
housing 

Mo and Cu 

150 (Mo) 
200 (Cu) 

14 Mo + 13 Cu 
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Vacuum chamber conditions were recorded periodically using ionization 
gauges located at the end of the chamber near the plane of the accel- 
erator grid and at the chamber wall si 2 m downstream from the thruster. 
The latter position is felt to be the most representative of the ambient 
vacuum chamber conditions. Figure 11 presents a time history of the 
vacuum chamber pressure recorded during a 14-hr erosion monitor test 

conducted at beam conditions of 1.5 A and 3000 V. After the beam 

“6 

voltage was applied, the tank pressure rose to the high 10 • Torr 
range as a result of ion bombardment of the beam collector. After 
this initial pressure rise, the chamber pressure dropped to a low point 
of 3.5 X 10~ Torr near the end of the 14-hr test. 

f . Post-Test Data Analysis 

Thruster operating times were scaled to remove several of 
the erosion monitor layers. After the tests, the monitors were removed 
from the various thruster components for analysis. To facilitate 
determining the number of layers removed, a "diagnostic etch" was 
performed on each monitor by ion-beam machining a groove through the 
remaining layers. By loosely positioning a mask above the surface, 
as shown in Figure 12, a tapered groove was machined to facilitate 
identifying layer interfaces at the boundary of the region masked 
during the test. 
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REGION MASKED 
DURING EROSION TEST 



DIAGNOSTIC ETCH REGION OF SPUTTERING EROSION MONITOR 

Figure 12. Diagnostic etch for erosion monitor analysis. 
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The erosion monitor layers, particularly those constructed with thin 
substrates, sometimes became blistered. Using thicker substrate 
material alleviated this problem somewhat by confining the blistering 
to regions near the edge of the baffle and at the tip of the cathode 
pole piece. This localization effect and the dependence on substrate 
thickness suggest that blistering may be caused by high temperatures. 
Blistering was observed only with monitors mounted on the baffle, 
baffle support, and cathode pole piece; none occurred with monitors 
mounted on the screen grid, accelerator grid, or neutralizer housing. 

Once the number of copper layers remaining had been determined, 
the erosion rate, W, was calculated using the following expression: 

W = [{N - N’ + 0.5) 5.ra,to ^ (N - N') «cu/VTa,Mo]/’' • 


where N and N' are the initial and final number of copper layers, 
respectively; 6 is the layer thickness; t is the test time; and y 
adjusts for the increased erosion rate of copper (y-p^ = 40 and 
^Mo ~ first tenn in the numerator represents the thickness 

of eroded Ta or Mo, with the factor 0.5 introduced to account for the 
unknown fraction of the Ta or Mo layer remaining above the highest 
visible copper layer. Adopting this definition of the unknown, 
remaining layer thickness, the uncertainty in the number of eroded 
layers is ±0.5. Thus the erosion rate uncertainty, U, can be expressed 
as 


^Ta,Mo - ^Ta,Mo/^'^ 


( 6 ) 


The second term in Eq. 5 accounts for the copper layers eroded, and the 
factor 1/y converts copper layer thickness to equivalent tantalum or 
molybdenum layer thickness. 

g. Results 

Table 5 presents thruster wear rates calculated from the 
multilayer erosion, monitor test results. The observed discharge 
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Table 5. Thruster Wear Rates Obtained froin Alternating-Film 
Erosion Monitor Measurements 


Parameters 


Case 2 

Case 3 

Case 4 

Operating Conditions 





Beam current, A 

2.0 

2.0 

1.5 

•1-Q- 

Beam voltage, kV 

1.6 

3.0 

3.0 

3.0 

Wear Rates, nm/hr 





Screen grid 





Center 

52.3 

52.3 

20.7 

2.6 

Mid- radius 

31 .3 

15.5 

9.1 

1.3 

Accelerator grid 





Center 

41.8 

36.5 

6.8 

3.9 

Mid-radius 

15.5 

15.5 

6.8 

1.3 

■ Baffle 





Discharge side 

76.8 

66.5 

33.8 

19.2 

Cathode side 

5.0 

5.5 

6.7 

3.8 

Baffle support 

■ 15.3 

15.3 

6.7 

3.8 ■ 

Cathode pole piece 

15.3 

15.3 

2.2 

1 .3 

Neutralizer 

167.5 

115.8 

107.9 

61 .3 

Uncertainty 

±5 

±5 

±3.2 

±1 .3 

^Neutralizer uncertainty 2.5 times the values shown. 
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chamber wear rates are relatively independent of beam voltage at the 
2- A beam current level. The screen grid wear at the center of the grid 
was unaffected by the variation in beam voltage, and the accelerator 
grid wear rate at the center was lower at the higher beam voltage. 

This effect is consistent with the increase in propellant utilization 
and the decrease in accelerator grid drain current that was observed at 
higher beam voltages. The wear rates are strongly dependent on beam 
current, which is to be expected because both single and double ion 
densities within the discharge chamber increase with beam current. 
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Multilayer erosion monitors attached to the downstream side of 
the accelerator grid indicated that the charge-exchange ion erosion 
pattern was asymmetric with respect to the accelerator grid hole 
•pattern. An indication of poor beamlet compensation was also evident 
in the angular dispersion data obtained from the E x B momentum 
analyzer. This is illustrated in Figure 13 » which shows the beamlet 
dispersion at various radial locations. The dispersion data at 
thruster radii of 0.5 and 0.75 indicates that the beamlets, were vectored 
off-axis by as much as 8 to 12®. 

F. THRUSTER PERFORMANCE AND SCREEN GRID EROSION MODELS 

This section presents mathematical models of thruster performance 
and screen grid erosion. The thruster performance model, which, is 
based on EMT performance documentation obtained previously at 3000 sec 
specific impulse in the 1 to 2.5 kW power range, gives an accurate 
mathematical description of the EMT performance parameters with beam 
current and beam voltage as independent variables. The model was used 
to calculate thruster performance parameters in the extended performance 
range (Figure 1); the results are compared with the experimental data 
obtained with the extended- performance thruster. An analytical model 
of screen grid erosion was derived, and erosion rates calculated with 
the model are compared with the wear rate data obtained using the 
multilayer erosion monitors. 

■ 1 . Performance Model 

The data from the extensive testing and performance documentation 
undergone by the 30-cm EMT was used to develop an analytical model of 
thruster performance. Standard curve fitting techniques were used to 
develop the mathematical description of EMT performance in terms of 
beam current, Ig, and beam voltage, Vg: 

Total power (in W) 

Pt = Ie(Vg -f- 200) + 65 , (7)- 
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CURRENT, ARBITRARY UNITS 



.Figure 13. Beamlet dispersion data obtained from E x B momentum 
analyzer measurements (r^ is the radius of the 
accelerator aperture array). 
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electrical efficiency 


Tl„ = 




e Pj ’ 


(8) 


propellant utilization efficiency (measured) 


^u = B 




(9) 


3 = 1 .0 - 0.8 


^BX 


2.2 + I 


BX J 


OO) 


'BX 


0 . % < 1 
^B “ ^ ’ ^B ~ ^ 


U A 


(ir) 


neutral loss rate (in A) 


I|^ = 0.240 + 0.032 Ig , 


( 12 ) 


thrust correction factor 


Y . 0.942 - 0.005 Ig + 

D 


(13) 


total thruster efficiency (corrected) 


TI4. = Y r\Ti 


e”u ’ 


(14) 
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specific impulse (corrected, in sec) 


^sp = 100 Y , (15) 

and thrust (corrected, in N) 

T - 2.039 X 10~^ Ylg ^ . (16) 

where 3 is the ratio of total propellant flow to the sum of the beam 
current and neutral loss rate, and is a dummy variable. The 
.expressions presented above were formulated to obtain a "best fit" to 
t'he EMT performance data in the 1- to 2-A beam current operating 
Ifjange. 

k The EMT performance model was used to calculate thruster performance 
parameters for the extended-performance range. These results, and the 
experimental results obtained during testing of the extended-performance 
thruster, are presented in Figures 14 through 18. The mathematical 
performance model gives an excellent description of the power, thrust, 
and jtotal efficiency of the extended-performance thruster. There is 
appreciable scatter in the experimental data for propellant utilization 
and consequently also for specific impulse. The data points correlate 
quite well with the analytically determined curves, although the 
points do not fall near the curves in all instances. Some of the 
discrepancy may have been caused by changes in beam voltage and magnetic 
baffle current since these quantities were constant in the EMT data 
used to formulate the model. 

Although beam current and beam voltage were treated as independent 
variables in the extended-performance operating range, this is true, 
only for the thruster operating envelope defined by (1) the perveance 
line of the ion optics, (2) the specified voltage standoff capability 
of the propellant isolators, and (3) the recommended maximum and 


41 



TOTAL POWER, kW 


6789-8 



0 0.5 1.0 1.5 2.0 2.5 

BEAM CURRENT, A 


Figure 14. Total power versus beam current: comparison of values 
predicted by the performance model with measured* 
experimental results. 
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THRUST, mN 
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TOTAL POWER, kW 


Figure 16. Total efficiency versus total power; comparison of 
values predicted by the performance model with 
measured experimental results. 
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Figure 17. Propellant utilization versus beam current: comparison of 

values predicted by the performance model with measured 
experimental results. 
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SPECIFIC IMPULSE, sec 



minimum beam current operating limits of 2.5 and 1.0 A. The perveance 
limit can be derived from the definition of perveance 

T 


’ (17) 

-5 -3/2 

where the perveance, p, is found from Figure 5 to be 2.4 x 10 AV . 

2. Screen Grid Erosion Model 

Thruster lifetime is a major concern in all electric propulsion 
applications because long thrust times are required. This is particularly 
true of the Halley's comet mission, in which thrusters may be required 
to operate for as long as 14,000 hr. The lifetime limitation of the 
current 30-cm EMT design is thought to be caused by failure of the 
screen grid due to sputter erosion damage. Based on the 10,000-hr 
endurance test results of thruster SN 701 , the erosion rate at the 
center of the screen grid was estimated at 35 nm/hr when operating at 
the 2 A beam current level. ^ At this rate, the time required to erode 
a substantial fraction (say 90%) of the screen grid is « 10,000 hr. 
Although this lifetime would be adequate for many missions, it is only 
marginal and perhaps inadequate for the primary propulsion require- 
ments of the Halley's comet rendezvous. A mathematical model of screen 
grid erosion is therefore needed to predict screen grid lifetime and 
to identify thruster modifications likely to reduce screen grid 
erosion. 

3. Sputtering Model 

The surface sputtering rate, W, can be expressed in terms of the 
incident ion flux, j, and the target sputtering yield, S, by 

W = (jV + 1 j+V^)/qtl , (18) 
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where n is the density of the target material, q is the' electronic 
charge, and the superscripts refer to the charge on the ion. This 
expression can be rewritten in terms of the measurable beam parameters 
by introducing the continuity equation for a nonreacting one-dimensional 
flow of the charged species 


J 



(19) 


where J is the current density immediately downstream from the accel- 
erator grid, and ratio of screen grid open area to 
active beam area. Combining Eqs. 18 and 19 yields an expression for 
the wear rate in terms of beam conditions: 


A 



1 + 


(20) 


where Jj is the sum of the single- and double-ion current densities. 
This expression can be applied if local values of the current density 
and sputtering yield (ion energy) are used. 

The screen grid lifetime can be predicted by applying the wear 
rate expression to the thruster centerline, where the total current 
density and double-to-single ion current density ratio are highest 
and maximum screen grid erosion is observed. The current density, Jjq, 
at the centerline is given by 



1 

F 


A. 


(21) 


where F is the beam flatness parameter, defined as the ratio of the 
average beam current density to the maximum beam current density 
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evaluated at the accelerator grid. Combining Eqs. 20 and 21 yields 
the screen grid erosion rate on centerline, 


W 


0 



qnF^B 


+ 0.5 ^ 
__ 



( 22 ) 


where is the ratio of double- ion current density to single-ion 
current density on centerline. Expressing the erosion rate in this 
form permits easy identification of those factors affecting screen 
grid lifetime.’ For example, operation at a high beam current increases 
the erosion rate and shortens thruster lifetime. This effect is non- 
linear because also increases with beam current. The inverse 
dependence of erosion rate on the product F<|>g indicates the advantages 
of a uniform current density profile and of a screen grid design having 
a high open area fraction. The advantages of thruster operation at 
reduced discharge voltage levels are implicit in Eq. 22 since the 
sputtering yields and ^ are strong functions of the discharge voltage.. 
This effect is also highly nonlinear since S^, and are nonlinear 
functions of discharge voltage. 


4. Thruster Lifetime Prediction 

The screen grid erosion expression can be written in a form con- 
venient for calculating thruster lifetimes by introducing the screen 
grid parameters and estimated beam flatness of the 30-cm EMT: 

= 0.75 

Ag = 573 cm^ 

22 -3 

n = 6.4 X 10^^ cm 


F = 0.58 . 
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Substituting these values into Eq. 22 yields 


Wq = 1.409 X 10 


S'^ + 0.5/ 
0 


1 + 


<^0 


Ig in nm/hr 


(23) 


+ ++ 

The sputtering yields S and S are evaluated at the single-ion energy 
E and the double-ion energy 2E using the experimental data of Askeroy 
and Sena. The single-ion energy E is assumed equal to the local 
plasma potential, which is approximated as the sum of, the discharge 
voltage plus the anode sheath potential (assumed to be| 6 V). Thus S 
can be approximated as 

E = Vp + 6 . '' (24); 


Eq.- 23 can be used to predict the thruster lifetime under various . . 
operating conditions, where the thruster lifetime, is defined as the' 
time required to erode a substantial fraction (90%) of the .screen grid. 
Table 6 presents calculated thruster lifetimes based on operating 
parameters taken from Table 2 and Ref. 3 and on the assumption of 90%. of 
the grid having been eroded. These results indicate the strong dependence 
of thruster lifetime on beam current operating point. Reducing -the beam 
current by a factor, of 0.5 results in an increase in predicted lifetime 
by a factor of almost 4. A factor of 2 is gained by the reducing beam 
current density; the remainder of the improvement occurs because the ratio 
of double-ions to single-ions is reduced. Table 6 also presents erosion 
rate and thruster lifetime calculations for discharge conditions that 
were measured using SHAG optics on the EMT thruster. At the 2-A beam cur- 
rent level, reducing the discharge voltage to 30 V by using SHAG optics 
increases the predicted thruster lifetime to 23,000 hr. this is a factor 
of 2.5 increase over the predicted lifetime when standard EMT grids are 
used. 

The validity of the thruster lifetime predictions can be assessed ’ 
by comparing the calculated screen grid erosion rates with measurements ' 
obtained during the endurance test of thruster SN 701. The measured 
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Table 6. Calculated Screen Grid Erosion Rates and 
Thruster Lifetimes 


Optics 

A 

\i 

^0 

X 10^ 

s"^ 

X 10^ 

“o- 

nm/hr 

3?, 

hr 

EMT 

2.0 

36 

0.5 

6.5 

0.4 

38 

9,000 

EMT 

1.5 

36 

0.3 

6.5 

0.4 

22 

15,300 

EMT 

1.0 

36 

0.15 

6.5 

0.4 

11 

31 ,500 

SHAG 

2.0 

30 

0.22 

3.5 

0.25 

15 

23,000 
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erosion at the center of the screen grid after 10,000 hr of operation 
at 1.4 A beam current was 0.18 mm, equivalent to an erosion rate of 
18 nm/hr. This value can be extrapolated (approximately) to 1.5 A beam 
current operation by multiplying by the ratio of beam currents. This 
results in an extrapolated screen grid erosion rate of 19.3 nm/hr. 

The screen grid erosion rate calculated for the extended-performance 
thruster operated at 1.5 A beam current is 22 nm/hr, which is about 
14% higher than the extrapolated value; this is considered to be 
reasonably good agreement.- The screen grid erosion rate of thruster 
SN 901 was 31 nm/hr afters; 4,000 hr of operation at 2 A beam current. 
The calculated value for the extended- performance thruster is 38 nm/hr. 
This agreement is considered reasonable, particularly considering the 
simplicity of the model and the uncertainties present. 

5. Comparison of Calculated and Measured Screen Grid Erosion Rates 

Calculated (using the erosion model) and measured (using the 
multilayer erosion monitor technique) screen grid erosion rates are 
presented in Table 7 for comparison. The agreement at the 1.5 A beam 
current level is within the uncertainty in the measured results and 
is therefore considered very good. Both the measured and calculated 
values are in agreement with the screen grid erosion rate of 19.3 nm/hr 
extrapolated (also for 1.5 A beam current) from the 10,000 hr endurance 
test of thruster SN 701. Although there is apparent agreement at the 
1.5 A beam current level, at higher or lower beam current levels the 
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Table 7. Comparison of Calculated and Measured 
Screen Grid Erosion Rates 


^B’ 

A 

Screen Grid Erosion Rate, nm/hr 

Calculated 

Measured 

1.0 

11 

2.6 ±1.3 

1.5 

22 

20.7 ±3.2 

2.0 

38 

52.3 ±5.0 
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agreement is poor. There are, however, several plausible explanations 
for the differences in calculated and measured erosion rates;' these 
are discussed below. 

The calculated results were based on sputtering yields taken from 
published data that had been taken in laboratory experiments conducted 
under carefully controlled vacuum conditions. Therefore, this sput- 
tering yield data may not be equivalent to what would have been obtained 
in the actual discharge chamber environment because of various effects 
such as surface temperature variations, adsorption of vacuum chamber 
residual gases, chemical reaction on discharge chamber surfaces, and 
formation of surface irregularities or "cones." Depending on the 
thruster power level and vacuum chamber conditions, these effects may 
either increase or decrease the effective surface sputtering yield. 
Another factor that strongly affects the calculated erosion rate is 
the beam flatness parameter. A value of F equal to 0.58 (taken from 
EMT beam measurements) was used in all calculations. However, the 
beam profile may become more uniform at reduced beam current levels, 
which would reduce the calculated screen grid erosion rate. The 
calculated erosion rate is strongly dependent on the shape of the 
beam profile: if the assumed value of beam flatness were reduced from 

0.58 to 0.47, the calculated erosion rate corresponding to the 2-A 
beam current conditions would be 47 nm/hr. This is within 10% of the 
measured value, which is less than the uncertainty in the erosion 
rate measurement. 
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The erosion rate measured at the 2-A beam current level is 
significantly higher than the 29 to 31 nm/hr measured during previous 
tests^’^ and, if valid, would severely limit the lifetime of the 
screen grid. One factor that may have increased the screen grid erosion 
rate (and, for that matter, the erosion rate for all discharge chamber 
components) is the low tank pressure maintained during the extended- 
performance thruster tests. The dependence of measured discharge chamber 
erosion rates on ambient conditions in the vacuum chamber was recognized 
early in the test program. This was based on an experimental and analyti- 
cal evaluation performed at NASA LeRC of the results from the 900 series 
thruster endurance testing at Hughes.^ Precautions were taken to ensure 
a high-quality vacuum environment during the erosion monitor tests. Other 
investigations of discharge chamber sputtering have indicated that baffle 
erosion is strongly dependent on cathode flow rate. In these studies, 
the baffle erosion rate was observed to increase by a factor of 2 when 
the cathode flow rate was increased from 80 to 135 mA. The cathode flow 
rate during the tests of the extended-performance thruster was about 
no mA at the 2-A beam level, which may explain the higher observed baffle 
erosion. Since the increase in baffle erosion with cathode flow is 
apparently caused by increased ionization near the baffle aperture, its 
effect should be localized. Therefore, it is difficult to associate an 
increase in screen grid erosion with cathode flow effects unless these ‘ 
flow effects give rise to some other phenomenon (such as increased dis- 
charge chamber oscillations) that strongly affects erosion rates. 

F. CONCLUSIONS 

Performance characteristics of the thruster used in this program 
are essentially identical to those previously documented and analytically 
modeled using the 30-cm EMT. Therefore, the mathematical model des- 
cribing the performance parameters of the EMT in terms of the independ- 
ent variables of beam current and beam voltage is applicable to the 
extended-performance thruster operating envelopes proposed for the 

Halley's comet rendezvous mission. 

s 
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Re-spacing the EMT ion optics assembly to an interelectrode gap 
of 0.76 mm provides more than adequate perveance throughout the 
extended performance operating envelopes and provides sufficient 
voltage standoff capability. The only other thruster modifications 
required for high-voltage operation are replacement of the cathode 
and main propellant isolators and the use of insulated restraints 
for preventing contact between wiring and thruster structural members. 

An accelerator grid voltage of 500 V is adequate to prevent 
electron backstreaming for beam voltages as high as 3900 V. The back- 
streaming limit of the re-spaced EMT ion optics is ~ 400 V. 

The ion beam envelope does not change significantly with beam 
conditions. When operating at a constant accelerator grid voltage of 
500 V, 95% of the primary beam ions are contained within a 20° cone 
and all primary beam ions are within a 51° cone. The ion beam cone 
angle can probably be reduced by as much as 10° by reducing the off- 
axis current component of those beamlets emanating from the outer 
grid apertures. 

High-power thruster operation results in component temperatures 
that are within the acceptable operating ranges. 

At 1.5 A beam current level, both the multilayer erosion monitor 
measurements and the screen grid erosion calculations agree with the 
screen grid erosion rate observed in the 10,000-hr endurance test of 
thruster SN 701. However, at higher or lower beam current levels, the 
agreement between calculated and measured results is poor. 

The measured discharge chamber erosion rates are substantially 
higher than those measured' in the EMT thruster; if these rates are 
valid, erosion would severely restrict thruster lifetime. The observed 
increase in measured erosion rates is thought to be due, at least in 
part, to carefully controlled vacuum conditions. 

An analytical model of screen grid erosion was developed that 
relates screen grid erosion to the. discharge chamber plasma properties 
and certain discharge chamber parameters. The erosion model suggests 
that the most effective means of improving screen grid lifetime are 
(1) to reduce discharge voltage by using small-hole accelerator grids 
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and (2) to improve the plasma density profile (increase the flatness 
parameter) near the screen grid by some means such as a magnetic field 
modification. 
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SECTION 3 

DEVELOPMENT OF THE HIGH-VOLTAGE PROPELLANT 
ELECTRICAL ISOLATOR 

The 900-series EMT was developed for operation at a specific 
impulse of 3,000 sec (beam voltage of 1,100 V). The propellant elec- 
trical isolators were developed to have a conservative, steady-state 
voltage rating of 1,500 V, As pointed out in Section 2, the extended- 
performance operation range includes beam voltages as high as 3,900 V. 
Consequently, a high-voltage propellant isolator design is required 
that can provide a conservative, steady-state voltage rating of at 
least 5,000 V. Because the allotted schedule time v/as short, the 
approach taken here was to scale existing isolator designs, build 
prototype components, and perform comparative tests. Two isolator 
design concepts were investigated: a multisection concept that is used 

Q 

in the 900-series EMT design and an insulating labyrinth concept that 
has- been successfully demonstrated in the British ion propulsion 
program. The evaluation of these concepts is described in this 
section. 

A. THE MULTISECTION HIGH-VOLTAGE ISOLATOR CONCEPT 

The fundamental problem in designing a propellant electrical 
isolator is the prevention of the ionization of the gaseous propellant 
being supplied to the thruster. Since the propellant vaporizers and 
propellant storage tanks should be operated at spacecraft potential, 
the full beam voltage appears across the insulator used to transmit 
the propellant to the thruster. Ionization of the gas in this situa- 
tion is governed by the type of gas, the gas pressure (or density), and 
the dimensions which characterize the vessel containing the gas. 
Empirical relations have been determined under controlled geometric 
configurations for several gases, including mercury vapor. These 
relationships, known as Paschen breakdown curves, relate breakdown 
voltage to the product of the gas pressure and the gap distance across 
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which the voltage is applied (p-d). A typical Paschen curve is shown in 
Figure 19. The basic physical principle involved- is that there exists 
a limiting voltage for any given product of gas pressure and charac- 
teristic dimension. Below that voltage, not enough electrons are 
accelerated across the dimension by the applied voltage to cause a 
cumulative ionization (avalanche); above that voltage, ionization 
occurs. This relationship between breakdown voltage and p-d is strongly 
dependent on boundary conditions (material, surface condition, etc.}. 
Measurements like those shown in Figure 19 are usually made using 
highly polished, plane, parallel electrodes to define the character- 
istic dimensions across which the voltage is applied (electric field is 
established). The voltage minimum on this curve in Figure 19 is 
sometimes referred to as the Paschen minimum. Voltages below this 
value cannot cause ionization of the gas in question at any value of p-d.- 
The EMT propellant electrical isolator design is based on the concept of 
a Paschen minimum breakdown voltage. 

In the 900-series EMT isolator design, the vapor flow channel is 
divided into seven chambers (or sections) by optically dense metal 
screens, as shown in Figure 20. The voltage is assumed to divide 
evenly across the length of the alumina cylinder that contains the 
vapor. Therefore, with 1500 V applied, the voltage across any section 
would be « 215 V, which is less than the Paschen minimum for mercury 
covered electrodes («250 V). At high flow rates, the breakdown voltage 
for this design is about 2000 V (see Figure 21). Because of the long 
lead time for obtaining high-quality ceramic materials, we elected to 
scale the- design conservatively (by a factor of four) to leave as many 
options as possible open. An isolator having 28 sections was therefore 
specified. One option is to maintain the electric field strength 
along the insulator surface constant and to increase the overall iso- 
lator length also by a factor of four. This would result in an assem- 
bly that is structurally questionable with regard to its ability for 
surviving vibration during launch. We elected to evaluate an assembly 
of this type (see Figure 22) to prove the scaling principle. The 
other design extreme is to keep the overall isolator length the same 
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Figure 19. 

Paschen curves for various gases. (Taken from Ref. 28.) 
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Figure 20. Isolator- vaporizer configuration 
of the 30-cm thruster. 
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Figure 21. Isolator breakdown voltage as a function of mercury flowrate. 
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Figure 22. 900-series isolator scaled directly 
to 28 sections (4 X). 
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while scaling only the internal elements. This approach is shown in 
Figure 23; the resulting component becomes a direct replacement for the 

isolator used on the 900-series EMT. The structural characteristics of 
this assembly should be the same as those of the 900-series design, but 
with the field strength along the insulator surface increased by about . 
a factor of four. Based on our experience with high-voltage insulators, 
surface breakdown on insulators is not significant until the field 
strength (in vacuum) exceeds 5 kV/cm. Similarly, the interfaces 
between the .dielectric and conductor require no special treatment (to 
avoid field concentrations) until field strengths exceed 5 kV/cm. 

Since the field strength in the design shown in Figure 23 is only about 
half this value, the design was considered to be adequately conserva- 
tive. The spacers selected for this design are thin alumina discs 
(such as shown in Figure 24). These discs, which are cast alumina, 
were chosen more for their short delivery schedule than for their 
material composition, although their insulating properties are essen- 
tially the same as those of the A1 300 alumina material used for the 
900-series isolator spacers. 

An isolator was constructed in a cathode-isolator-vaporizer 
configuration for each of-these designs; the same quality-control 
procedures were used as are used in fabricating isolators for the 
30-cm EMT. Both isolators were successfully tested with .applied 
voltages up to 6 kV without breakdown under full vapor flow conditions. 
Test procedures and detailed results are described in Section 3.C. 

B. INSULATING LABYRINTH ISOLATOR CONCEPT 

With this concept, the propellant vapor is forced to flow through 
an insulating labyrinth. This can be done by using (1) a porous 
ceramic that has relatively small pores, (2) tightly packed insulating 
spheres, or (3) tightly packed insulating fibers. The physical prin- 
ciple is to prevent ionization by providing a high surface- to-volume 
ratio so that any electron- ion pair formed will reach a surface before 
it collides with another gas atom. This approach was used successfully 
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^ALUMINA SPACERS 0.65 mm (0.025 in.) THICK (SEE FIGURE 20). 


Figure 23. 


High-voltage propellant electrical 
section design concept. (Internal 


isolator using multi - 
seal ing.) 
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2.880 cm (1.134 in.) diameter 
x2.883 cm (1.135 in.) diameter 


0.32 cm (0.125 in.) 
diameter holes 
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0.063 ± 0.005 cm 
(0.025 + 0.002 in.) 


Figure 24. Alumina spacer for the 28-section isolator. 
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in the British thruster development programJ^ Although the design 
criteria for this isolator concept are not well defined, our approach 
was to use an isolator envelope that had the same dimensions as the 
EMT isolator and to replace the spacers and screens by tightly packed 
ceramic spheres, as shown in Figure 25. We assumed that we could use 
Rye's results, which show breakdown voltage to be inversely proportional 
to sphere diameter. Given that relationship, spheres of 0.1 to 
1 mm diameter should have been used. But the only insulating spheres 
that were obtainable during this program were a mixture ranging from 
0.2 to 2 mm in diameter. We elected to assemble an isolator with the 
configuration shown in Figure 25 in a manner such that the downstream 
boundary of the isolator could be removed for replacement of the 
insulating spheres. With such an isolator, we proposed to select sphere 
sizes from the assortment and evaluate the dependence of breakdown 
voltage on sphere diameter. 

A demountable main isolator assembly (as shown in Figure 25) was 
successfully fabricated. But the sample of ceramic spheres obtained 
was found to contain sphere sizes that were predominantly in the 
1.2 to 1.5 mm diameter range. Consequently, the evaluation performed 
is neither complete, nor conclusive. The test procedures and results 
are discussed in more detail in Section 3.C. 

C. ISOLATOR TESTS 

After fabrication, the isolators were inspected visually, given 
a preliminary electrical checkout for heater continuity, and installed 
in a vacuum bell jar. The vacuum facility is pumped by a 15-cm- 
diameter mercury diffusion pump and has both freon refrigerated and 
liquid nitrogen baffles. Typical ambient pressure is about 2 x 10 ® Torr. 
A schematic of the test apparatus is shown in Figure 26. Isolator test- 
ing proceeded as follows: 

• Isolator leakage current was 'measured with a voltage of 
6 kV applied; the isolator and vaporizer were at ambient 
temperature. 
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Figure 25. High-voltage propellant electrical isolator using insulat- 
ing labyrinth design concept. 
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VAPORIZER AUXILIARY 

HEATER POWER HEATER 




ISOLATOR 
HEATER POWER 


Figure 26. Schematic of high-voltage propellant isolator 
test apparatus. 
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• Isolator leakage current was measured with 6 kV applied; 
the isolator base was heated to 300°C. 

• With the isolator base at 300°C, the vaporizer flow rate 
was adjusted to wl50 mA (equivalent); voltage was applied 
to the isolator in 1-kV increments and leakage current 
was measured as a function of applied voltage up to a 
limit of 6 kV or to isolator breakdown, whichever occurred 
first. 

The results of this testing are summarized in Table 8. Both of the 
multi section isolators meet the operational temperature, voltage, and 
propellant flow requirements (300°C, 6 kV, and 100-to 150-mA propellant 
flow rate); the insulating-labyrinth-type isolator sustained internal 
breakdown as a function of isolator temperature. Since the short inter- 
nally scaled multi section isolator appeared to operate as well as did the 
longer version, it was’ selected for a longer duration test of 300 hr. 

The conditions for the 300-hr test were specified by NASA: 

• Isolator base temperature = 300°C 

• Vaporizer temperature = 375°C 

• Applied voltage = 4 kV. 

The isolator leakage current was monitored continuously by strip chart 
recorder and flow rate data was recorded twice each working day. The 
average flow rate measured for the test was about 260 mA (approximately 
three times the EMT flow rate). The leakage current rose approximately 
linearly with time over the duration of the test, as shown in Figure 27. 
Leakage current behavior similar to this has also been observed in 
testing performed during the development of the 900-series EMT isolator 
design. The leakage is related to surface contamination of the 
ceramic insulator. This contamination causes a semiconducting thin 
film to form on the insulator; this film. conducts in proportion to 
the absolute temperature of the isolator. (A detailed description of 
both the sources of contamination and preventive measures is given in 
Ref. 12.) In brief, the insulator surface must be protected from any 
possible contamination during fabrication by careful assembly procedures 
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Table 8. Summary of Isolator Test Results 


Conditions for Isolator Leakage 
Tests 

Long 

riulti section 
Isolator^ 

Short 

Multi section 
Isolator® 

Isolator 
Filled with 
Ceramic 
Spheres'^ 

Ambfent temperature; 6 kV applied 
Isolator leakage, nA 

2 

3 


3 

Isolator hot; 6 k\/ applied 
Isolator base temperature, °C 

290 

300 ' i 

'290 

Vaporizer temperature, 

80 

190 

260 

Isolator leakage, yA 

0.13 

0.02 j 

0.1 

06 

Isolator hot, vaporizer 
flow rate "normal" 

Isolator base temperature, °C 

330 

330 

302 

152 

Vaporizer temperature, °C 

345 

340 

255 

258 

Vaporizer flow rate, A 

0.150 

0.130 

1.81 

2.0 

Isolator leakage at 1 kV, yA 

0.09 

0.02 

0.005 

0.1 

Isolator leakage at 2 kV, yA 

0.1 

0.03 

0.01 

0.2 

Isolator leakage at 3 kV, yA 

0.1 

0.03 

c 

0.3 

Isolator leakage at 4 kV, yA 

0.11 

0.03 


0.4 

Isolator leakage at 5 kV, yA 

0.12 

0.04 


0.54 

Isolator leakage at 6 kV, yA 

0.13 

0.04-ld 


e 


^These isolators were cathode-isolator-vapori 2 er configurations 
■operated with a cathode attached for flow impedance. 

K * 

This isolator was a main- isolator- vaporizer configuration . 
^Breaks down at 2800 V. 

^Meter reading fluctuated between 0.04-1 yA — no evidence of 
breakdovi/n. 

^Breaks down at 5200 V. 
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Figure 27. • Isolator leakage current versus elapsed test time for 
the multisection isolator design shown in Figure 23. 
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and during operation by eliminating possible contaminants from the 
isolator vicinity and by including protective shielding. Procedures 
prescribed by NASA LeRC have eliminated this leakage for EMT operating 
temperatures and voltages. Although the prescribed handling and 
operating precautions were observed in preparing this test, portions 
of the auxiliary heater (which was used to simulate thruster thermal 
input) operated at too high a temperature (700°C) and were too close 
to the isolator. This could have been a source of contamination. 

Another possible source of contamination is sputter deposition on 
the insulator from micro-discharges between the protective shields. 
Although discharges could not be seen visually, the electrometer was 
observed to fluctuate sporadically at applied voltages above about 
3.5 kV. Although determining the source of the observed leakage current 
was beyond the scope of this study, we recommend further investigation 
to determine whether a problem actually exists. 

D. CONCLUSIONS 

The study described in this section was a "single-pass" evaluation 
of propellant electrical isolator technology. Although refinement of 
the designs evaluated could be expected to affect the results quantita- 
tively, we feel that the following general conclusions can be drawn: 

• The multisection isolator approach can be scaled for higher 
voltage operation based on the Paschen minimum voltage 
principle. 

• An isolator with adequate voltage capability can be 
demonstrated in a configuration that does not perturb 
the design of the EMT. 

• The insulating labyrinth isolator approach can provide 
high-voltage isolation, but the design criteria are not 
well established. 

The problem of isolator leakage at the high voltages used (4 to 6 kV) 
is unresolved, but should be tractable by the same procedures used in 
developing the EMT isolator. Therefore, we conclude that the design 
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and development of a high-voltage propellant electrical isolator to 
meet the requirements of the Halley's comet mission can be accomplished 
within the required time period and at a Tow technical risk. 
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SECTION 4 


A CHARGE EXCHANGE PLASMA EFFLUX MODEL 

One PMaC approach studied under this program would provide the 
high voltage necessary to operate the ion thrusters by a direct series 
connection of the solar cells in the solar array. Consequently, large 
areas of the solar panel must operate at high positive potential with 
respect to the rest of the spacecraft. Although the high-energy ion 
beam exhaust is not expected to be influenced by these high-potential 
surfaces, the low-energy ions that will form near the thruster by 
charge exchange could present a problem. These ions drift out of the 
beam, driven by thermal velocities and potential gradients in and 
around the edge of the high-energy ion-beam plume. The electrons 
may possibly be conducted from the thruster neutralizer through the 
dilute "plasma" formed by these charge-exchange ions and be collected 
by the high-voltage solar array. Our objective was to develop an 
analytic model to describe the charge-exchange ion efflux from a 30-cm 
EMT. This model was developed and presented to NASA LeRC during the 
first month of the program for use in analyzing possible interactions 
between the spacecraft and the solar array that could result from 
operating the ion thrusters using a high-voltage solar array. 

The charge-exchange plasma model presented here was derived pri- 
marily on the basis of experience and measurements made with a single 
30-cm thruster. Simplifying assumptions were made to permit analytic 
descriptions to replace empirical data. The results obtained are 
considered to be as accurate as any published in the literature. No 
attempt was made to include assymmetries associated with operating a 
single neutralizer. Similarly, it was considered to be an overspeci- 
ficationto infer the variation in space potential outside of the ion 
beam. The charge exchange ion production was calculated and it was 
assumed that a reasonable potential distribution existed within the 
ion beam. This produces an ion density and velocity distribution at 
the beam edge (the beam edge is defined here as the boundary of the 
envelope containing 95^ of the primary ion beam current). 
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A. CHARGE-EXCHANGE ION CALCULATIONS 

Several simplifying assumptions were used in calculating the 
charge exchange ion flux: 

• The neutral atom efflux is emitted in a cosine distribution 
with a uniform distribution at z = 0. The total neutral 
atom efflux was 250 mA (equivalent current if each atom 
were charged). 

• The distribution of the primary beam current is Gaussian 
with a spreading half angle of 150 . A beam current of 

2 A was used. 

• The charge exchange ions travel in a radial direction and 
are emitted from a conical boundary enclosing 95% of the 
primary beam current. 

• The potential distribution within the ion beam region is 
described by the Boltzmann equation with an electron 
temperature of 5 eV. 

The neutral and primary ion beam density distributions are described 
below in more detail, followed by the charge-exchange calculation. 

1 . Neutral Density Analysis 

Figure 28 shows the geometry used to calculate the downstream 

neutral density. In these calculations it was assumed that the 

neutrals are emitted uniformly over the thruster diameter at z = 0 

with a cosine distribution. This assumption permits calculating the 

neutral efflux passing through an annular region having radii r-j and 

ro a distance z downstream from the thruster by using formulas used 
2 11 .. 
for calculating radiant heat transfer. The view factor F-j 2 is 

defined as the fraction of the total flux leaving surface 1 (assumed 

to emit in a cosine distribution) that lands on surface 2. The neutral 

efflux Hg^ leaving the thruster and landing on the annular region is 

thus 





(25) 
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z =0 


SURFACE 2 



Figure 28. Geometry used for neutral density calculations. 



2 

Since the flux is equal to nvA, where n is the density (parti cl es/cm ), 
V is the velocity (cm/sec), and A is the area (cm ), the neutral 
density at the thruster, n^, is related to the density at the annular 
region n (r-j,r2,z) by 


n^ A-| ^12 ^2 


(26) 


or 




(r^ 


7^ ^12 


where 


^12 p ^ 2 " ’"l^ ^ 

^ *"o 

^2 = ^ '"2 * - ‘”'2 >'0 

T, = , ,2 ^2 , 


(27) 


(28) 


(29) 


(30) 


As r^ ^ r^, Eq. 27 represents the neutral density n(r,z) at radius r 
where r = 1/2 [r^ + r^). This relationship is plotted in Figure 29 as 
a function of r with Z as a parameter. These curves show how the 
square distribution at the thruster becomes progressively flatter as 
z increases. Numerically, for a total neutral efflux of 250 mA (equiva- 

2 4 

lent), an area of 660 cm , and a neutral velocity of 2.4 x 10 cm/sec 
('x. 300°C) we have 


0.250 

660 x 2.4 X 10^ X 1.6 X 10 


-19 


10 


11 atoms 


(31) 


cm 
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This value together with the distribution shown in Figure 29 and 
described by Eq. 27 were used for all the charge-exchange calculations. 
Since the charge-exchange generation rates are directly proportional 
to n^^, the generation rate at any other neutral efflux value can be 
scaled proportionally. This neutral efflux for the 30-cm EMT is 
essentially independent of ion beam current in the range of currents 
considered {1 to 2.5 A). 

2. Primary Beam Current Density 

At any z position, the primary current density was assumed to 
have a Gaussian profile with a radial dependence described by 

'JCr) . (32) 

Since the total beam current is given by 


!■*■=/ 2TrrJ(r,z) dr = vrl T , (33) 

0 ■ 

the current density on axis (r = o) is 

Jq " 

The current contained up to radius r is 

r 

I(r) = / 2irrJ(r,z) dr (35) 

0 


or 
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2 


I(r) = I"" 



(36) 


L J 

When r equals ,(r) = [l - e ^] = 0.95 (i.e., 95% of the 

beam current is contained in a radius defined as r^c = r„)- This 

yb 0 

rgg radius is defined to be the beam "edge" at which the charge-exchange 
current and density was calculated. In, reality, there is no beam edge 
where the beam density abruptly drops to zero. It was- further assumed 
that the rnn current envelope spreads with a 15° half angle. This 

4 * + • 

beam divergence is representative for I = 2 A, V = 1100 V, and 
V" = 500 V. From the above assumptions it follows that 


rgg = 15 + z,sin 15° 


cm 


(37) 


The resulting current-density profiles calculated using Eqs.'25, 27, 

and 30 are shown in Figure 30 where J(r,z) is plotted versus r with 

z as a parameter. The ordinate in Figure 30 can be converted to density 

in ions/cm^ by multiplying by 1.92 x 10^ (i.e., 1 mA/cm^ at 1100 V 

is 1.92 X 10® ions/cm^). For a 1 A beam at 3.9 kV, multiply by 

0.5 y 11 00/3900 X (1.92 x 10®) = 5.07 x 10®; for a 2.5 A beam at 

2.75 kV multiply by (2.75/200) x 71100/2750 x (1,92 x 10®) = 1.67 x 10®. 

3. Charge Exchange Calculations 

The charge exchange current can be easily calculated from the 
expressions for the primary ion current J(r,z) and neutral density 
n(r,z). The total charge exchange current dl created in a path 
length dz is 

CO 

^ [j(r,z) 2irr. n(r,z) dr] a dz , (38) 

CX Q 
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where a is the charge exchange cross section, which was taken as 
7 X 10“^^ cm^ from Ref. 14. This is its value at 1100 V; it increases 
slowly with decreasing voltage. The total current created in the 
space between the thruster and the axial location, 0 to z, is given by 

= / «cx • <39) 

This expression was numerically integrated, and the result is shown 
in Figure 31. Most of the charge-exchange current is created close to 
the thruster; for example, 75% of the total charge-exchange current is 
created within one thruster diameter and 90% is created within two 
diameters. These calculations of the total generated current do not 
depend on the position of the beam edge or on the direction in which the 
ions move. 

B. CHARGE EXCHANGE ION EFFLUX 

Although the calculation of the charge-exchange ion-generation 
rate does not depend on the position of the beam edge or on the direc- 
tion in which the ions move, the calculation of the ion efflux from 
the beam does. The high-energy ions are essentially mono-energetic 
and have velocities that are directed primarily along the thrust 
axis. If there were no electrons injected into this beam, a logarithm- 
ically varying potential gradient would be established in a radial 
direction (the beam edge would expand rapidly with increasing distance 
from the plane of the accelerator). -Ideally, the neutralizer supplies 
electrons to equal the ion space charge density and thus there is no 
net charge at any point in the neutralized ion beam envelope. But in 
reality, the electrons have thermal energy and their velocities are 
randomly directed and higher than the velocities of the highly directed 
heavy ions. Consequently, a radially directed potential distribution 
is established to prevent electron loss at the beam edge (as determined 
by the ion trajectories from the outermost electrode apertures). In 
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AXIAL POSITION DOWNSTREAM FROM THE THRUSTER, cm 


Figure 31. Total charge-exchange ion formation in interval from 

z = 0 to distance ^ downstream of thruster. (Note that 
75% of the total charge-exchange ions generated are 
formed within a distance of one thruster diameter 
(30 cm)}. 
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where a is the charge exchange cross section, which was taken as 
7 X 10 cm^ from Ref. 14. This is its value at 1100 V; it Increases 
slowly with decreasing voltage. The total current created in the 
space between the thruster and the axial location, 0 to z, is given by 

This expression was numerically integrated, and the result is shown 
in Figure 31. Most of the charge-exchange current is created close to 
the thruster; for example, 'v 75% of the total charge-exchange current is 
created within one thruster diameter and 90% is created within two 
diameters. These calculations of the total generated current do not 
depend on the position of the beam edge or on the direction in which the 
ions move. 

B. CHARGE EXCHANGE ION EFFLUX 

Although the calculation of the charge-exchange ion-generation 
rate does not depend on the position of the beam edge or on the direc- 
tion in which the ions move, the calculation of the ion efflux from 
the beam does. The high-energy ions are essentially mono-energetic 
and have velocities that are directed primarily along the thrust 
axis. If there were no electrons injected into this beam, a logarithm- 
ically varying potential gradient would be established in a radial 
direction (the beam edge would expand rapidly with increasing distance 
from the plane of the accelerator). Ideally, the neutralizer supplies 
electrons to equal the ion space charge density and thus there is no 
net charge at any point in the neutralized ion beam envelope. But in 
reality, the electrons have thermal energy and their velocities are 
randomly directed and higher than the velocities of the highly directed 
heavy ions. Consequently, a radially directed potential distribution 
is established to prevent electron loss at the beam edge (as determined 
by the ion trajectories from the outermost electrode apertures). In 
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AXIAL POSITION DOWNSTREAM FROM THE THRUSTER, cm 


Figure 31. Total charge-exchange ion formation in interval from 

z = 0 to distance z dov/nstream of thruster. (Note that 
0 . 75% of the total charge-exchange ions generated are 
formed within a distance of one thruster diameter 
(30 cm)). 
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the absence of physical boundaries, electron motion along the beam 
axis is unrestricted and any significant potential gradient (relative 
to the radial potential gradient) would be rapidly eliminated by 
electron motion. But if electrons are able to escape from the beam 
in an .axial direction, then an axial potential gradient will also be 
established. This would be the case for an ion thruster operated in 
a vacuum chamber. Consequently, experimentally measured axial varia- 
tions in potential are expected to depend strongly on the boundary 
conditions in the experiment, and we do not consider the existence of 
an axial potential gradient to be representative of operating an ion 
thruster in space. Therefore, to calculate the charge-exchange current 
density we will assume that the ions move in a radial direction and that 
the charge-exchange current density differential dJ^^(r,z) is defined by 




dj (n,z) = - 9 — - 
cx 2 tt r, 


95 


This is a worst case assumption since a small amount of current 
generated outside r^g has been included. This current density is 
plotted in Figure 32 and is shown in tabular form in Table 9. Fig- 
ure 32 shows the rapid decrease in 0 with distance. For example, at 
one thruster diameter downstream the current density is 12 % of the 
value at z = 0. To calculate the charge-exchange density at rg^ 
requires estimating the average velocity in the radial direction at 
the beam edge, which is a function of the beam plasma potential dis- 
tribution, V(r). 

Kaufman^ ^ has shown that a good approximation to the potential 
distribution of the beam is described by the Boltzmann equation that 
relates the plasma density, potential, and electron temperature as 
follows: 


eV 
" kT, 

n = no e 


(40) 
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Table 9. Tabular Data for Charge-Exchange Data 
Plotted in Figures 5 and 6 


z 

(cm) 

rg5 

(cm) 

•^cx, 

[iA/cm^ 

n 

ions/cm^ 

102 

41.40 

1.04E-01 

2.07E+07 

104 

41.91 

9 . 86E-02 

1.97E+07 

106 

42.43 

9.39E-02 

1.88E+07 

108 

42.95 

8.96E-02 

1.79E+07 

110 

43.47 

8.54E-02 

1.71E+07 

112 

43.98 

8.16E-02 

1.63E+07 

114 

44.50 

7.79E-02 

1.56E+07 

116 

45.02 

7.45E-02 

1.49E+07 

118 

• 45.54 

7.13E-02 

1.43E+07 

120 

46.05 

6.82E-02 

1.36E+07 

122 

46.57 

6.54E-02 

1.31E+07 

124 

47.09 

6.26E-02 

1.25E+07 

126 

47.61 

6.01E-02 

1.20E+07 

128 

48.12 

5.77E-02 

1.15E+07 

130 

48.64 

5.54E-02 

l.llE+07 

132 

49.16 

5.32E-02 

1.06E+07 

134 

49.68 

5.11E-02 

1.02E+07 

136 

50.20 

4.92E-02 

9.83E+06 

138 

50.71 

4.73E-02 

9.46E+06 

140 

51.23 

4.55E-02 

9.11E+06 

142 

51.75 

4.39E-02 

8.77E+06 

144 

52.27 

4.23E-02 

8.45E+06 

146 

52.78 

4.07E-02 

8.15E+06 

148 

53.30 

3.93E-02 

7.86E+06 

150 

53.82 

3.79E-02 

7.58E+06 

152 

54.34 

3 . 66E-02- 

7.32E+06 

154 

54.85 

3.53E-02 

7.07E+06 

156 

55.37 

3.41E-02 

■ 6.83E+06 

158 

55.89 

3.30E-02 

6 . 60E+06 , 

160 

56.41 

3.19E-02 

6.38E+06 

162 

56.92 

3.08E-02 

6.17E+06 

164 

57.44 

2.98E-02 

5.97E+06 

166 

57.96 

2.89E-02 

5.78E+06 

168 

58.48 

2.80E-02 

5.59E+06 

170 

58.99 

2.71E-02 

5.42E+06 

172 

59.51 

2.62E-02 

5.25E+06 

174 

60.03 

2.54E-02 

5.09E+06 

176 

60.55 

2.47E-02 

4.93E+06 

178 

61.06 

2.39E-02 

4.78E+06 

180 

61.58 

2.32E-02 

4 . 64E+06 

182 

62.10 

2.25E-02 

4.50E+06 

184 

62.62 

2.19E-02 

4.37E+06 

186 

63.14 

2.12E-02 

4.24E+06 

188 

63.65 

2.06E-02 

■ 4 . 12E+06 

190 

64.17 

2.00E-02 

4.00E+06 

192 

64.69 

1.95E-02 

3.89E+06 

194 

65.21 

1.89E-02 

3.78E+06 

196 

65.72 

1.84E-02 

3.68E+06 

198 

66.24 

1.79E-02 

3.58E+06 


Z 

(cm) 

^95 

(cm) 

•Icxl 

yA/cm2 

n 

ions/cm3 

2 

15.52 

1.29E+01^ 

2.59E+09 

4 

16.03 

1.13E+01 

2.26E+09 

6 

16.55 

9 . 66E+00 

1.93E+09 

8 

17.07 

8.18E+00 

1.64E+09 

10 

17.59 

6.91E+00 

1.38E+09 

12 

18.10 

5.84E+00 

1.17E+09 

14 

18.62 

4.95E+00 

9.91E+08 

16 

19.14 

4.22E+00 

8.44E+08 

18 

19.66 

3.61E+00 

7.23E+08 

20 

20.18 

3.11E+00 

6.22E+08 

22 

20.69 

2.69E+00 

5.39E+08 

24 

21.21 

2.34E+00 

4.69E+08 

26 

21.73 

1.05E+00 

4.10E+08 

28 

22.26 

1.80E+00 

2.60E+08 

30 

22.76 

1.59E+00 

3.18E+08 

32 

23.28 

1.41E+00 

2.82E+08 

34 

23.80 

1.26E+00 

2.52E+08 

36 

24.32 

1.12E+00 

2.25E+08 

38 

24.83 

l.OlE+00 

2 . 02E+08 

40 

25.35 

9.10E-01 

1.82E+08 

42 

25.87 

8.23E-01 

1.65E+08 

44 

26.39 

7.46E-01 

1.49E+08 

46 

26-90 

6.79E-01 

1.36E+08 

48 

27.42 

6.19E-01 

1.24E+08 

50 

27.94 

5.67E-01 

1.13E+08 

52 

28.46 

5.20E-01 

1.04E+08 

54 

28.97 

4.78E-01 

9.55E+07 

56 

29.49 

4.40E-01 

8.80E+07 

58 

30.01 

4.06E-01 

8.13E+07 

60 

30.53 

3.76E-91 

7.52E+07 

62 

31.04 

3.49E-01 

6.97E+07 

64 

31.56 

3.24E-01 

6.47E+07 

66 

32.08 

3.01E-01 

6.02E+07 

68 

32.60 

2.81E-01 

5.61E+07 

70 

33.11 

2.62E-01 

5.24E+07 

72 

33.63 

2.45E-01 

4.90E+07 

74 

34.15 

2.29E-01 

4.59E+07 

76 

34.67 

2.15E-01 

4.30E+07 

78 

35.19 

2.02E-01 

4.04E+07 

80 

35.70 

1.90E-01 

3.80E+07 

82 

36.22 

1.79E-01 

3.58E+07 

84 

36.74 

1. 68E-01 

3.37E+07 

86 

37.26 

1.59E-01 

3.18E+07 

88 

37.77 

1.50E-01 

3.00E+07 

90 

38.29 

1.42E-01 

2.84E+07 

92 

38.81 

1.34E-01 

2.69E+07 

94 

39.33 

1.27E-01 

2.55E+07 

96 

39.84 

1.21E-01 

2.42E+07 

98 

40.36 

1.15E-01 

2.29E+07 

100 

40.88 

1.09E-01 

2.18E+07 


^Expressed in scientific notation, e.g. , 1.29E+01 = 1.29 x 10^, 
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Figure 32. Calculated current density of Hg charge-exchange ions 
passing through conical surface enclosing 95% of beam 
current as a function of distance downstream from the 
thruster. 
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For this equation to be consistent with 


it is necessary that the following equation hold: 


eV(r) = kT 


At rnc where r = \/T r , V(r) = 3 kT^ = -15 V, where the center of the 
95 ^0 e 

beam at r = 0 is at V = 0. Let g(r,z) be the volume charge-exchange 
generation rate (rin A/cm }, where g(r,z) equals J(r,z) n{r,z) cr. Then 
the average velocity <u>, at the rgg beam edge is given by 


<u(z)> = g(r,z) - V(r)]- 2irr g.{r,z) 2irr dr .'(43) 


This expression was evaluated at several z positions from z = 0 to 
z = 200 cm. The average energy is essentially independent of z with 
a value 10.1 ± 2 eV (which is 2 kT^). The corresponding average 
velocity, <v>, can therefore be obtained by solving 


1 m<v>^ = 2 kT^ 

2 e 


which yields 


<v> = 2.\ ^ = 3.1 X 10^ cm/sec 
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which is twice the Bohtn minimum velocity for formation of a stable 
sheath (the Bohm velocity was derived for a monoenergetic velocity 
distribution). Using this value of <v>, the charge-exchange current 
densities were converted to ion number densities using n = J/e<v> and 
plotted in Figure 33 as a' function of distance downstream from the 
thruster, 

C. EFFLUX MODEL FOR A THRUSTER ARRAY 

To estimate the charge-exchange boundary for an array of thrusters, 
the current-density distribution for a 6 x 2 thruster array was calcu- 
lated at z = 30 cm and z = 60 cm downstream from the array. The pri- 
mary ion beam current density distribution at an arbitrary point was 
obtained by summing the current density contributions for each of the 
12 thrusters at the observation point r^: 


where 1 is the distance between the center of the thruster 

and the observation point. Contours were calculated for normalized' 

-1 -2 -3 

amplitudes of e , e , and e . The latter is the 5% amplitude 
defined to be the nominal beam edge that encloses 95% of the beam 
current. The results of these calculations are shown in Figure 34, 

The contours shown in Figure 34 show that at z = 30 cm the thrust 
beams have not merged (i.e., the source boundary for the charge 
exchange ions are circular) and that at z = 60 cm the rg^ boundaries 
have merged to form an undulating rectangular boundary. For either 
case, however, the far-field current density distribution must be 
obtained by calculating trajectories that pass through the beams of 
adjacent thrusters; such a calculation would be beyond the scope of 
this work. 
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n (r, z), CHARGE-EXCHANGE ION DENSITY, lONS/cm 
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Figure 33. Calculated density of Hg charge-exchange ions at conical 
surface enclosing 95% of beam current as a function of 
distance downstream from the thruster. 
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A rather crude estimate' can be made for combinations of thrusters 
by assuming that the charge-exchange space charge will collect in the 
region between thrusters until a potential distribution is formed such 
that the neutralizing electrons are confined. At this point, charge 
exchange ions will again be driven out of the beam region at the same 
rate at which they are generated. Figure 35 shows two examples of com- 
bined beams at z = 60 cm. The current densities and shown in 
Figure 35(a) refer to the curved and linear segments of the combined 
beam envelope. The current density, , is simply the value given in 
fable 9 for z = 60 cm. 

The current density ^2 given by 


J 


2 


60.5fr X J-j 

2 X 53.3 = '^1 • 


(47) 


For the 12-thruster case (Figure 35(b)), J-j and J2 are obtained in 
the same manner; however, is given by 


_ 10 X 61 X IT , 

3 2 X 320 1 



(48) 


Hence the geometry and current density for a combined beam envelope can 
be obtained using the above expressions and Table 9 for any value of z. 


D. THE HOLLOW CATHODE NEUTRALIZER 

% 

The hollow cathode neutralizer used on the 30-cm ion thruster 
is fundamentally an electron source that provides electrons to the ion 
beam to equalize or neutralize the ion current. The configuration in 
use is shown as Figure 36 with densities and potentials as indicated; 
The relatively low-voltage (12 to 20 V) discharge that is maintained 
between the cathode orifice and the keeper (keep alive) electrode 
serves to maintain a high-density plasma in the cathode orifice region 
and to supply electrons to the ion beam. If the neutral and charged 
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(b) TWELVE THRUSTERS OPERATING ■ 


Figure 35. Assumed cross section for combined ion beams at z 
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Figure 36. Neutral izer/i on-beam diagram. 
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particle densities in the orifice and surrounding regions are as shown 
in Figure 36 , the process of beam neutralization is thought to be 
relatively effective and any azimuthal assymmetries that might be 
present are very localized to a region within several cm of the neutral- 
izer. A complete description of hollow cathode and neutralizer opera- 
tion can be found in Ref. 16 through 27. 

E. CONCLUSIONS 

The generation of charge-exchange ions in the ion-beam exhaust 
plume takes place predominantly in the region immediately downstream 
from the accelerator electrode. The neutralized ion beam establishes 
a radial potential gradient to contain the neutralizing electrons and 
this potential gradient expels the charge exchange ions from the ion 
beam in a radially outward direction. A beam envelope boundary was 
defined and the charge-exchange ion density and ion current density were 
calculated as a function of location on the boundary. The results 
considered to be valid provided that the strength of externally applied 
electric fields are significantly smaller at the beam boundary than the 
field required to contain the thermalized electrons (from the neutral- 
izer) in the beam. We presented a technique for scaling the charge- 
exchange ion efflux in proportion to the ion beam current and also 
described a procedure for using the calculated data from a single 
thruster to obtain the efflux from a thruster array. This data, 
together with the geometry assumed, is considered to form a representa- 
tive numerical model for the charge-exchange ion efflux that is useful 
for estimating the plasma densities in the space around the thrust 
system. The ion beam and efflux model described should also be useful 
for formulating the numerical boundary conditions necessary for com- 
puting the electron and ion currents that will flow between the ion 
beam plume and any spacecraft or solar array surfaces that would be 
operated at high relative potentials (as in the case of the direct 
drive solar array). 
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SECTION 5 
CONCLUSIONS 


Each section of this volume contains in detail the conclusions 
reached for the particular technology area’ investigated. Therefore, 
here we only summarize the important conclusions: 

• The 900-series 30-cm EMT requires only minor modifications 
to satisfy the requirements of a Halley's comet rendezvous 
or other extended performance missions. 

• The performance characteristics of the thruster that was 
modified for extended-performance operation were the same 
as had been projected from extrapolated EMT data. 

• A high-voltage propellant electric isolator design that is 
based on the EMT isolator concept can be scaled to meet 
the extended-performance requirements (as defined under 
this study). 

® The capability of EMT technology (the technology includes the 
modifications required for extended-performance operation) 
for meeting the lifetime requirements of an extended 
performance mission was not demonstrated and requires 
verification. 
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